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Abstract: We complete the physical picture for the evolution of a high-energy jet propagating 
through a weakly-coupled quark-gluon plasma by investigating the thermalization of the soft com¬ 
ponents of the jet. We argue that the following scenario should hold: the leading particle emits a 
signihcant number of mini-jets which promptly evolve via quasi-democratic branchings and thus de¬ 
grade into a myriad of soft gluons, with energies of the order of the medium temperature T. Via elastic 
collisions with the medium constituents, these soft gluons relax to local thermal equilibrium with the 
plasma over a time scale which is considerably shorter than the typical lifetime of the mini-jet. The 
thermalized gluons form a tail which lags behind the hard components of the jet. We support this 
scenario, hrst, via parametric arguments and, next, by studying a simplified kinetic equation, which 
describes the jet dynamics in longitudinal phase-space. We solve the kinetic equation using both 
(semi-)analytical and numerical methods. In particular, we obtain the first exact, analytic, solutions 
to the ultrarelativistic Fokker-Planck equation in one-dimensional phase-space. Our results confirm 
the physical picture aforementioned and demonstrate the quenching of the jet via multiple branching 
followed by the thermalization of the soft gluons in the cascades. 
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1 Introduction 

It is by now well established that, within weakly-coupled QCD at least, the energy loss by an energetic 
parton and/or the associated jet propagating through a dense medium, such as a quark-gluon plasma, is 
dominated by medium-induced radiation, that is, the additional radiation triggered by the interactions 
between the partons from the jet and the medium constituents [1-11] (see also the review papers [12- 
14]). This picture has led to a rather successful phenomenology, based on, or at least inspired by, 
calculations in perturbative QCD, which has allowed one to understand many interesting observables 
at RHIC and the LHC, like the nuclear modification factor or the suppression of di-hadron azimuthal 
correlations in ultrarelativistic nucleus-nucleus collisions [15-18]. 

More recently, the experimental studies of the phenomenon known as ‘di-jet asymmetry’ in Pb-|-Pb 
collisions at the LHC have demonstrated that a substantial fraction of the energy loss by an energetic 
jet is carried by relatively soft hadrons propagating at large angles with respect to the jet axis [19-26]. 
This pattern too can be understood, at least qualitatively, within the pQCD picture for medium- 
induced radiation, which predicts the formation of well-developed gluon cascades, or ‘mini-jets’, via 
multiple branching [27-35] (see also Refs. [36-39] for earlier, related, studies and the recent review 
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paper [40]). Within these cascades, the energy is efficiently transmitted, via quasi-democratic branch¬ 
ings, from the ‘leading particle’ — the parton that has initiated the (mini)jet — to a large number 
of comparatively soft gluons, which can be easily deviated towards large angles by rescattering in the 
medium. The theoretical description of medium-induced multiple branching started being developed 
only recently and in its current formulation it leaves unanswered a number of important questions. 

The basic question is, what is the microscopic mechanism responsible for the energy loss ? Of 
course, one already knows that, from the perspective of the leading particle, the dominant mechanism 
at work is radiative energy loss and that the energy carried by the primary radiation is efficiently 
transmitted to softer and softer particles, via successive branchings in the gluon cascades. But what 
is the physical mechanism which stops these cascades, and at which energy scale (i.e., what is the 
low-energy end of the cascade) ? How does this mechanism influence the dynamics of the branchings 
at higher energy scales ? And where does the energy go, when it flows out of the cascade ? 

In order to better appreciate these questions, it is useful to briefly recall our current understanding 
of the medium-induced jet evolution via multiple branching (see Sect. 2 below for details). The gluon 
cascades which are relevant for us here are those generated by iterating gluon emissions of the BDMPSZ 
type^ [2-5]. The BDMPSZ mechanism governs the emission of relatively hard gluons, with energies 
to ^ T, which undergo multiple scattering in the surrounding medium. Here, T is the characteristic 
energy scale of the medium, say, the temperature for the case of a quark-gluon plasma in thermal 
equilibrium, or, more generally, the average pT of the background hadrons. The distinguished feature 
of this mechanism for our present purposes is the fact that it favors ‘quasi-democratic branchings’, 
that is, 1 —)• 2 gluon splittings where the daughter gluons carry comparable fractions of the energy of 
their parent gluon. Such branchings occur fast and are extremely efficient in redistributing the energy 
among the branching products: they lead to turbulent cascades, in which the energy flows from one 
parton generation to the next one, without accumulating at intermediate steps. 

In most theoretical analyses so far, on has assumed, for simplicity, that the branching dynamics 
remains unmodified down to arbitrarily low energies. This allowed for elegant and physically trans¬ 
parent solutions [28, 30, 35], which exhibit wave turbulence with a characteristic scaling spectrum — 
the analog of the Kolmogorov-Zakharov spectrum [41, 42] for the medium-induced cascades — and 
have interesting consequences for the energy loss by the jet (see Sect. 2 below). But such an ‘ideal’ 
cascade leads also to unphysical results: after a finite interval of time (the ‘branching time’ t]^j.(E), to 
be specified in Sect. 2) the original energy E of the leading particle gets transmitted to quanta which 
are arbitrarily soft {uj —)• 0) and hence can propagate at arbitrarily large angles. While this peculiar 
‘final state’ is clearly unacceptable, it is important to stress that the phenomenon of wave turbulence 
is in fact more general and could very well coexist with a physically acceptable final state: the scaling 
spectrum survives unchanged if one stops the branching process at some finite energy scale p* <C ill, by 
introducing there a ‘perfect sink’. The ‘perfect sink’ — a concept familiar in the theory of turbulence 
[41, 42] — is, by definition, a mechanism which is capable to absorb the energy flux generated by the 
cascade at p*, without influencing the branching dynamics at higher energies w 3> p*. 

On physical grounds, there is an obvious candidate for such a ‘sink’ : the surrounding medium. 
The soft gluons from the cascade with a; ~ T are expected to thermalize via collisions in the medium 
and thus deposit their energy inside the medium. This motivated proposals in the literature to 
terminate the cascade at the medium scale T, e.g. by enforcing an ‘infrared’ cutoff p* ~ T on the 

^The acronym ‘BDMPSZ’ stands for Baier, Dokshitzer, Mueller, Peigne, Schiff, and Zakharov. 
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branching dynamics [30, 32]. But such previous arguments were insufficiently developed; it was not 
clear, e.g., what is the actual thermalization mechanism, why should this inhibit the branching process, 
and whether this should act as a ‘perfect sink’, or, on the contrary, wash out the wave turbulence. 
[Numerical simulations using an ad-hoc infrared cutoff p* [32] observed a strong distortion of the 
scaling spectrum, due to the accumulation of gluons in the bins above p*. With increasing time, this 
pile-up extends up to high energies cu » p* (see also the discussion in Sect. 5.2 below).] 

It is our main purpose in this paper to clarify these and related questions, via a dedicated theo¬ 
retical analysis. Specifically, assuming the medium to be a weakly-coupled quark-gluon plasma with 
temperature T, we shall study the possibility that the soft components of the jet, with energies a; ~ T, 
thermalize via elastic constituents with the quarks and gluons from the plasma. To that aim, we shall 
consider a special kinetic equation, which emerges via specific approximations from more general (but 
also more difficult to solve) equations existing in the literature [36, 43], and which will be argued to 
capture the interesting dynamics to parametric accuracy at least. This equation, to be introduced in 
Sect. 3, describes the evolution of the gluon distribution f{t, z, pz) created by the jet in the longitudinal 
phase-space, with the z axis referring to the direction of propagation of the leading particle. 

The kinetic equation includes two types of collision terms; an inelastic one, describing multiple 
branching with the BDMPSZ splitting rate, and an elastic one, which describes 2 —)> 2 collisions with 
the medium constituents in the Fokker-Planck approximation [44]. The latter is a priori suitable for a 
probe particle, or a dilute system of such particles, which can be distinguished from the thermal bath 
(like a heavy quark [45, 46]), but can also be applied to the gluons from the jet with relatively large 
momenta pz — oj T [47]. Indeed, the occupation numbers for the gluons generated via multiple 
branching remain very small down to a; = T, as we shall see (cf. the discussion in Sect. 2.1). 

We shall further argue that, as a result of thermalization, the branching process effectively termi¬ 
nates at the medium scale a; ~ T and we shall mimic this by inserting an ‘infrared’ cutoff p* ~ T in the 
splitting rate. The physical mechanism beyond this cutoff will be clarified too — this is related to the 
non-linear effects in the total gluon distribution, i.e. the distribution produced by the medium plus the 
jet (see the discussion in Sect. 3) — , but a proper treatment of this mechanism would require working 
with the non-linear kinetic equation obeyed by the full distribution, a task which is extremely hard in 
practice. As we shall see in our numerical simulations, the physics around this cutoff is smeared by 
elastic collisions and thermalization, so in practice we do not expect strong artifacts related to p*. 

Finally, the restriction to the longitudinal dynamics is needed to simplify the problem (in partic¬ 
ular, in view of numerical calculations) but it is also physically motivated. The longitudinal momenta 
remain much larger than the respective transverse components {pz p±) so long as a; ^ T, that is, 
during most stages of the dynamics, where they control the relevant time scales. This approximation 
fails, strictly speaking, in the approach towards thermal equilibrium, but even in that case it captures 
the correct time dependence to parametric accuracy. (In fact, a similar approximation has been used 
in all the previous studies of the in-medium cascade, including those which have explicitly considered 
the dependence upon transverse momenta [29, 31, 32, 34].) Furthermore, the gluon distribution along 
the longitudinal axis z is the most interesting one in view of a study of thermalization, since this is 
strongly inhomogeneous to start with: in the absence of collisions, all the gluons in the jet, even the 
softest ones, would propagate along the light-cone at z = t, together with the leading particle. 

In view of the above approximations, we expect that a thermalized distribution emerging from our 
kinetic equation should look like a ‘tail’ lying well behind the front of the jet, i.e. at \z\ <C t, where it is 
quasi-homogeneous, and which in longitudinal momentum features the classical thermal distribution 
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for massless particles in one spatial dimension, that is, the Maxwell-Boltzmann distribution 
This is indeed what we shall observe in our solutions. 

Specifically, we shall study the kinetic equation via a combination of (semi) analytic and numerical 
methods. For the Fokker-Planck dynamics alone, we will be able to find exact, analytic, solutions — 
in particular, the solution corresponding to a steady source along the light-cone (in Sect. 4.1) and the 
exact Green’s function in longitudinal phase-space (see Sect. 4.2.1). To our knowledge, these are the 
first examples of exact solutions for a relativistic Fokker-Planck equation. By combining this Green’s 
function with known, analytic, solutions for the (ideal) branching process [28, 30, 35], we have been 
able to give a relatively simple, semi-analytic, analysis of the full dynamics, under the assumption 
that the medium acts as a perfect sink at the scale p* = T. (Under this assumption, it is indeed well 
justified to treat the branching part of the dynamics as a source which injects gluons at the scale p* 
at a rate known from the previous analyses of the turbulent cascade; see the discussion in Sect. 4.) 

In order to test the ‘perfect sink’ assumption and also to have more a complete study which 
includes the interplay between branchings and elastic collisions at energies u > p*, we present in 
Sect. 5 a detailed numerical analysis of the kinetic equation, with infrared cutoff p* = T in the 
branching integral. The numerical solutions turn out to be qualitatively similar to the semi-analytic 
ones in Sect. 4, but they bring additional clarifications, in particular, on the role of the infrared cutoff 
p* and, related to that, on the limitations of the ‘perfect sink’ approximation. 

The overall physical scenario which emerges from these explicit solutions is in agreement with the 
general picture anticipated in Sect. 2, but it is more precise and also more complete than the latter. 
It can be summarized as follows (see Sects. 2 and 5 for details): 

A ‘leading particle’ (LP) with initial energy E ^ T which crosses the medium along a dis¬ 
tance L radiates abundantly (i.e., with a probability of order one) relatively soft gluons with energies 
w < Uhr{L). Here, Wbr(A) ~ o^qL?', with q ~ ln(l/as) the jet quenching parameter, is the 

characteristic energy scale for democratic branchings: gluons with smaller energies u < ti;br(A) can 
undergo a democratic branching within a time At < L, where those with higher energies uj ^ a;br(T) 
cannot. In the experimental conditions at the LHC, one typically has T <C a;br(L) E, so the LP 
belongs to the second category above, whereas the ‘primary gluons’ emitted by it belong to the first 
one. Accordingly, each of these ‘primary gluons’ generates a gluon cascade (‘mini-jet’) via successive 
democratic branchings. Each such a cascade ends at the thermal scale T, meaning that the energy 
w of a primary gluon gets distributed among a large number lo/T ^ 1 of soft gluons. These gluons 
undergo elastic collisions with the medium constituents and thus relax to a thermal distribution in 
momentum after a time of order tj-ei ~ T'^/q- At the same time they separate in z from the harder 
partons (which keep propagating along the light-cone) and thus form a tail at \z\ < t, which lags 
behind the front. Therefore the momentum distribution looks very different near the front of the jet 
(z ~ t), where the spectrum shows an approximate scaling behavior, as expected for ‘ideal’ branching, 
and in the tail at z < t, where the distribution is nearly thermal. 

The energy carried away by the gluons in the thermalized tail is naturally interpreted as the 
energy lost by the jet to the medium. This is controlled by the hardest ‘mini-jets’, those with energies 
oj ~ a;br(L), and therefore it scales with the medium size like L?. More precisely, the L^-scaling 
is the respective prediction of the ideal branching process, where the energy loss is computed as 
the energy carried by the turbulent flow [28, 30]. For the full dynamics including elastic collisions 
and thermalization, we numerically find that the energy deposited by the jet inside the medium is 
somewhat smaller than this ‘ideal’ prediction, albeit comparable to it. This reduction reflects the fact 
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Figure 1. The phase-space energy density \p\f/T produced by an energetic jet with E = 90T at two successive 
times: (a) an early time t = 0.1tbr(^')) when the jet is almost unquenched; (b) a larger time t — 0.3tbr(-E’)) 
when the jet is partially quenched. The secondary peak visible around p ~ T in plot (b) represents the energy 
lost towards the medium via thermalization. The reference scale tbr(i') is the characteristic time scale for the 
evolution of the jet via democratic branching, to be explained in Sect. 2.1 (this is given by Eq. (2.5)). 


that the medium is not a ‘perfect sink’, more precisely the fact that the soft gluons cannot thermalize 
instantaneously (see the discussion in Sect. 5.3.2). 

The above considerations are illustrated by the plots in Fig. 1 which are in fact extracted from our 
numerical results in Sect. 5. These plots show the (normalized) energy distribution {\p\/T)f{t,z,p) 
(with p = pz) produced by an incoming leading particle with E = 90T at two successive times, an 
early time, when the jet is almost unquenched, and a later time, when the jet is partially quenched. 
In the ‘late time’ plot in Fig. 1 (b), one clearly sees the accumulation of particles at soft momenta 
p < T, due to branchings and elastic collisions. 

Let us finally stress that the above picture strictly holds for a very energetic jet, with initial 
energy E ^ a;br(L), which is only ‘slightly quenched’ — meaning that the LP survives in the final 
state and the energy lost towards the medium is only a small fraction of the initial energy E. In that 
case, the final distribution, as visible in Fig. 1 (b), may be viewed as the superposition of the LP peak 
together with the distributions separately created by all the ‘mini-jets’. But the individual ‘mini-jets’ 
have lower energies uj < a;br(L), so they are fully quenched by the medium — their whole energy 
gets transmitted to the thermalized tail and the mini-jets disappear in the medium. The distribution 
created by a single mini-jet will be discussed too in Sect. 5. 

2 The physical picture 

In this section, we summarize the physical picture underlying the in-medium evolution of a jet gener¬ 
ated by a high-energy parton propagating through a weakly-coupled quark-gluon plasma. This picture 
largely reflects the current understanding of this problem, as emerging from the literature, but includes 
some additional arguments which will be physically motivated, together with some expectations to 
be subsequently confirmed by the new analysis in this work. We start with a brief review of recent 
studies of the medium-induced gluon cascade [27-30], which recognized the importance of multiple 
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branchings, but did not address the important problem of the thermalization of the soft components of 
the jet. Then we discuss the interplay between multiple branching and the elastic collisions responsible 
for thermalization. 

2.1 Multiple branching and the medium-induced gluon cascade 

An energetic ‘probe’ parton (gluon or quark) which propagates through a dense QCD medium, such 
as a quark-gluon plasma, undergoes elastic collisions with the constituents of the medium, leading to 
‘collisional’ energy loss — on the average, the energy transferred from the probe to the medium is larger 
than in the opposite direction — and to the broadening of the probe distribution in (longitudinal and 
transverse) momentum, due to the random nature of the ‘kicks’. Besides, the collisions trigger gluon 
emissions by the probe, leading to additional, ‘radiative’, energy loss, which in practice dominates 
over the collisional one, in spite of the fact that the emission probability is suppressed by a factor of 
as- This is possible because the coherence effects inherent in a quantum emission lead to a stronger 
dependence of the energy loss upon the medium size — the radiative component rises, roughly, like 
(with L the distance traveled by the probe through the medium), whereas the collisional component 
rises only like L. Notwithstanding, the elastic collisions will play an important role for the subsequent 
discussion — as we shall see, it is responsible for the in-medium deposition of the energy lost via 
radiation in a typical event. 

The precise mechanism responsible for medium-induced radiation depends upon the ratio between 
the gluon ‘formation time’ tform ~ (the typical duration of the emission, as fixed by the uncer¬ 

tainty principle) and the mean free path Amfp ~ l/[asTln(l/Q!s)] between two successive small-angle 
collisions. Here, uj is the energy of the emitted gluon and p± is its transverse momentum, as acquired 
during the formation time, via collisions. For a single collision, one typically has p\ ~ ~ ctsT’^-, 

with m£) the Debye mass. For a series of independent collisions occurring during a time interval 
At ::?> Amfp, one has p\ ~ ^At, where q ~ m^/Amfp ~ ln(l/Q;s) is the ‘jet quenching parameter’ 
— a transport coefficient which characterizes momentum diffusion. Using these estimates, one finds 
that, with increasing cj, one interpolates between a single-scattering (or ‘Bethe-Heitler’) regime at low 
energies uj < T, where tform ~ ujjmj^ < Amfp, and a multiple-scattering (or ‘LPM’, from Landau, 
Pomeranchuk, and Migdal) regime at high energies cu >> T, where tform ~ S> Amfp- In the LPM 

regime, there is a large number of collisions which coherently contribute to a single emission. 

In what follows, we shall consider the LPM regime alone. Indeed, we shall later argue that 
the branching process terminates around w ~ T, hence the phase-space for Bethe-Heitler radiation 
is comparatively small. The calculation of the gluon branching rate in the LPM regime to leading 
order in pQCD has been first given by Baier, Dokshitzer, Mueller, Peigne, and Schiff [2, 3, 6], and 
independently by Zakharov [4, 5]. One has thus obtained the following result for the differential 
probability per unit time and per unit x for the collinear splitting^ of a gluon with energy uj into two 
daughter gluons with energy fractions x and 1 — x, with 0 < x < 1: 

*I ^br _ ^ /2 

dxdt 27r ^£ 0 ^( 3 ^,^-^) 

^The splitting is effectively collinear since the transverse momentum squared acquired by the gluon during the for¬ 
mation time tform is much smaller than the respective quantity acquired after formation, namely during the lifetime 
tbr ~ tform/d of the gluou Until its next splitting (see below). 
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In this equation, Pg^g{x) = Nc[l — x(l — x)]^/x(l — x), with Nc the number of colors, is the leading 
order gluon-gluon splitting function of the DGLAP equation and t{orm{x,u) is a more precise estimate 
for the formation time, which involves the average energy x(l — x)uj of the two daughter gluons: 


^form(3:, Co) — 


I x(l — x)a; 


qes{x) = q[l-x{l- x)] 


( 2 . 2 ) 


qesix) 

Using Eq. (2.1), one can evaluate the probability for a branching to occur during a given interval At: 


AV 


~ X 


, d%, 

dx dt 


At 


~ OL . 


x(l — x')<jJ 


At, 


(2.3) 


where a = UsNc/tt. This probability becomes of order one, meaning that multiple branching is 
important during At, provided 


x(l — x)a; ~ a^qAt^ = a;br(At). 


(2.4) 


This condition can be satisfied by two types of emissions: 

(a) very asymmetric splittings, for which either x<Clorl—x<Cl, whereas the energy a; of the 
parent gluon can be relatively hard (for instance, oj S> a;br(At)); 

(b) ‘quasi-democratic’ branchings, where the two daughter gluons carry comparable fractions of 
the total energy, x~l — x~0(l), but the parent gluon is relatively soft: oj rsj Wbr(At). 

Reversing the argument for case (b) above, we conclude that it takes a time At ~ tbr(^^)) with 



for a gluon with energy u to undergo a ‘quasi-democratic’ branching [28, 36, 48]. This duration tbr(^^) 
should be compared to the medium size L which is available to that parton: 

(i) If tbr(w) L, then the parton with energy oj can emit abundantly — i.e. with probability of 
0(1) — only relatively soft gluons with x <C 1, in such a way that xu < ujhr(L) = a^qlP'. Accordingly, 
the original parton will ‘survive the medium’ : it will be recognizable in the final state due to the fact 
that its final energy will be considerably higher than for all the other gluons, as produced by radiation. 

(ii) If thr{^) ^ L, then the parton with energy cu will ‘disappear inside the medium ’ — it will 
undergo a quasi-democratic branching before it exits the medium and thus it will be replaced by a 
pair of softer gluons, which can ‘democratically’ split again. Eventually, the original parton will leave 
behind it a gluon cascade generated via successive, quasi-democratic splittings. Note that each new 
gluon generation in that cascade has a lower energy and hence a shorter lifetime than the previous 
ones. Accordingly, the overall lifetime of the cascade is of the order of the branching time tbr(w) of 
the initial gluon. 

Case (i) is the interesting situation for the leading particle (LP) which initiates a typical jet 
measured in Pb-I-Pb collisions at the LHC. Indeed, denoting the energy of this LP by E, one generally 
has E > 100 CeV, whereas the characteristic energy for multiple branching is much smaller: a;br(L) = 
o^qL^ ~ 12 CeV for a medium with q = ICeV^/fm and L = 5fm (we used d = 0.3). The above 
estimate also shows that, for the interesting values of L, the branching scale oj\y^{L) is much harder 
than the medium temperature T ~ 0.5 CeV; hence, one typically has T <C uj\,riL) ^ E. 
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Case (ii) applies to the typical primary gluons — the gluons which are directly radiated by the 
LP in a typical event —, which have relatively soft energies uj < and thus have the time to 

develop gluon cascades (‘mini-jets’) via quasi-democratic branchings. Harder primary emissions, with 
energies up to^ cUc = are possible as well (provided E > uic, of course), but these are rare events 
which occur with a probability of 0{a) and do not generate mini-jets. Such hard but rare emissions 
will not be explicitly considered in what follows, since they do not contribute to the energy loss by 
the jet towards the medium. (But they are important for the average energy loss by the LP [2-6].) 

The quasi-democratic nature of the splittings has important consequences for the energy flow 
across a mini-jet: it leads to wave turbulence [28, 30]. Via successive branchings, the energy flows 
from one gluon generation to the next one, without accumulating at any intermediate value of oj. The 
precise mathematical condition for wave turbulence is that the energy flux — the rate for energy flow 
along the cascade — should be independent of u. This condition is indeed satisfied for the gluon 
cascade generated by a primary gluon with energy wq < a;br(L), at least at sufficiently small values 
w <C Wo [28, 30]. (See the discussion in Sect. 4.2.2 for more details.) 

If this branching dynamics was to remain unmodified down to arbitrarily small values of w, then 
the whole energy would end up into a ‘condensate’ at w = 0 [28]. In reality though, we expect the 
gluon cascade to terminate at the thermal scale T, because the very soft quanta with w < T can 
efficiently thermalize via elastic collisions (see the discussion in the next subsection). If the medium 
acts as a perfect sink at the lower end of the cascade — in the sense of absorbing all the quanta with 
to < T without modifying the dynamics of branching at higher energies w S> T — then the whole 
energy carried by the turbulent flow is eventually transmitted to the medium. Under this assumption, 
the total energy loss by the jet towards the medium is obtained as [28, 30] 

ALlflow - ^Wbr(L) = , (2.6) 


with V ~ 4.96. This result, which is independent of the initial energy E of the LP and grows with 
the medium size like L^, truly applies (under the ‘perfect sink’ assumption) for a very energetic jet 
with E ^ a;br(L). In the opposite limit where E < a;br(L) (the case of a mini-jet), one rather has 
ALlfiow — E (see Eq. (4.30) for a more general expression). Eq. (2.6) admits a natural physical 
interpretation [30]: the LP radiates an average number u/2 of primary gluons with energies of order 
Wbr(L), which then transmit their whole energy to the medium, via democratic branchings followed by 
the thermalization of the soft gluons (w ~ T) at the lower end of the cascades. We shall later discover 
that the assumption that the medium acts as a ‘perfect sink’ has some limitations in practice, yet it 
can be used for qualitative considerations and parametric estimates. 

Another important assumption that was implicitly postulated by previous analyses of multiple 
branching [28, 30] is that the branching dynamics is linear, meaning that the gluons from the cascade 
can split, but not also recombine with each other. This assumption is correct provided the partons 
cascade are sufficiently dilute. The precise condition is that f{t,x,p) <C 1, where f{t,x,p) is the 
gluon phase-space occupation number (below, Ng the total number of gluons in the jet), 


f{t,x,p) 


(27r)3 dNg 
2(A(? — 1) (E’xd^p 


(2.7) 


®This upper limit cUc on the energy of medium-induced gluon emissions follows from the condition that the gluon 
formation time fform(ti^) be at most as large as L. 
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This quantity has not been explicitly computed in the previous studies, but it is straightforward to 
obtain an order-of-magnitude estimate for it via physical considerations. 

Consider a typical mini-jet, as generated by a primary gluon with initial energy wq < 

Over a time interval of order tbr(‘^o)) this whole energy gets redistributed, via multiple branching, 
among a large number Ng ~ wq/T of soft quanta with energies u) ^ T. Their occupancy can therefore 
be estimated as 

^ 

where pz ^ oj ^ T, Ap^ is the transverse momentum squared acquired by a gluon via rescattering in 
the medium, /S.z is the longitudinal extent of the distribution, and Ax^ is the corresponding spread 
in the transverse plane. 

The transverse phase-space Ap^^Ax^ occupied by a gluon with energy a; ~ turns out to be 
independent of a;. Indeed, during a time At, a gluon accumulates a transverse momentum broadening 
/S.p\ ~ gAt, leading to an uncertainty 


2 A 2 

Ax I ~ —^Ar ~ ^-TT- 


pI 




(2.9) 


in its transverse location. Taking At of the order of the gluon lifetime. At ~ tbr(w), one finds 


Ap^ Ax5_ — 




( 2 . 10 ) 


This is independent of a;, as anticipated, and parametrically large. Eq. (2.10) is truly a lower limit, 
since the soft gluons can also inherit part of the transverse momentum of their harder parents. The 
above argument also shows that, so long as co ^ T — as is the case for the gluons which control 
the splitting process — the transverse momenta remain much smaller than the longitudinal ones: 


P±<^Pz^ w. 

Consider now the longitudinal distribution of the soft quanta within the mini-jet. As already 
discussed, soft gluons are emitted promptly, so they can be produced anywhere along the cascade. 
After being emitted, they efficiently lose energy and randomize their direction of motion, via elastic 
collisions. Accordingly, they separate from each other and also from the harder (p ^ T) partons in 
the mini-jet over a time interval ~ treh which is small compared to the overall lifetime tbr(wo) of 
the cascade. Hence, on the average, these soft gluons should be quasi-uniformly distributed along z, 
within a distance /S.z ~ tbr(wo). Using this, together with (2.8) and (2.10), we finally deduce 
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This number increases with wq, so it is interesting to evaluate it for the largest possible value ujq ~ 
^hr{L) — the one which also controls the energy loss by the overall jet (cf. Eq. (2.6)). In that case, 
ibr(wo) ~ L, so one finds 


f{T) ~ 


a® qL 

Af r2 ■ 


( 2 . 12 ) 


The corresponding estimate for the jet as a whole can be simply obtained by multiplying this result by 
a number ~ v, cf. Eq. (2.6). The occupation number (2.12) is parametrically small at weak coupling 







and furthermore suppressed at large Nc- It is easy to check that the condition /(T) <C 1 is always 
very well satisfied in practice. 

The estimate in Eq. (2.12) should be more properly viewed as a lower limit on /: as we shall 
explain in Sect. 4.2.2, the soft gluons with w ~ T are more abundantly produced during the late stages 
of the cascade, at times t ~ so their distribution in z is not really homogeneous. (This will 

be also confirmed by the numerical simulations in Sect. 5; see in particular Fig. 8.) An upper limit on 
the longitudinal occupancy is however easily obtained by assuming the smallest possible longitudinal 
spread for the soft gluons, namely their lifetime tbr(T). With Az ~ tbr(T) and ojq ~ a;br(T), Eq. (2.8) 
implies (recall that tbr(T) ~ l/d^T), 


, . ^br(A) qLP' 

^ ^ THbr{T) ^ iV2 T 


(2.13) 


which is still much smaller than one in all practical situations of interest, as one can easily check. 

For what follows, one should keep in mind that the previous estimates for / refer exclusively to the 
soft gluons generated by the jet via multiple branching. When cu ~ T, these gluons add to those from 
the background medium, whose occupation numbers are given by the usual, Bose-Einstein, thermal 
distribution and hence are of order one. So the present results also show that the effect of the jet on 
the occupancy of gluons with oj <T represents only a small perturbation. This in particular implies 
that the jet cannot produce ‘hot spots’ in the medium (it cannot significantly increase the local energy 
density): the energy density carried by the soft jet constituents with w ~ T is obtained by multiplying 
/(T) by N‘^T^ and hence remains much smaller than the respective quantity for the thermal gluons, 
i.e. e ~ N^T^, so long as f{T) <C 1. 


2.2 Elastic collisions and thermalization 

In the discussion in the previous subsection, we have implicitly assumed that the only effect of the 
elastic collisions between the gluons from the jet and the medium constituents is to trigger new 
branchings. In reality though, such interactions can also transfer energy and momentum between 
the colliding particles, leading to energy loss and momentum broadening for the jet constituents. As 
we shall see, these effects remain negligible so long as the gluons in the cascade have relatively large 
energies uj ^ T, but they become a leading-order effect, and the driving force towards thermalization, 
when OJ T. 

To propose a theoretical description for these interactions, it is essential to recall, from the previous 
discussion, that the gluon system produced via multiple branching is dilute. Hence, it is appropriate to 
study the effects of elastic collisions on individual gluons from the jet. So long as the gluon under study 
has a relatively large energy w ^ T, it can be unambiguously distinguished from the thermal gluons, 
so it is possible to describe its dynamics by using the same methods as for other energetic probes, like 
a relativistic heavy quark (see e.g. [44-47]). In this subsection we shall use a Langevin description 
because of its formal simplicity. Later on, we shall employ the equivalent method of the Fokker-Planck 
equation for more elaborate studies. Strictly speaking, both descriptions will eventually fail when the 
energy of the ‘probe’ gluon decreases down to a; ~ T, but even in that case they remain qualitatively 
correct, in that they capture the correct time scale for thermalization to parametric accuracy. 

For more clarity, let us first assume that the branching dynamics is switched off, meaning that the 
energetic gluon suffers only elastic collisions in the plasma. The Langevin equation which encompasses 
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the effects of these collisions reads as follows: 


^ r, {e{t)eit')) = 1 -1'), (2.14) 

where u* = p^'/p, with i = 1,2,3, is the particle velocity, p is a friction coefficient, and is a 
stochastic force (the ‘noise’). Microscopically, the total force in the r.h.s. of Eq. (2.14) represents 
the Lorentz force generated by random, quasi-classical, fields in the plasma — the color fields of the 
thermal particles, which are slightly disturbed out-of-equilibrium by their scattering with the external 
particle. The ‘drag force’ /* = —r/v® describes the average effect of this microscopic force, which is 
the energy transfer from the probe to the medium, whereas the noise term represents its random 
component leading to momentum broadening. The average over the noise reflects the thermal average 
over the microscopic sources of this force. The facts that the noise correlator is local in time (‘white 
noise’) and also isotropic are non-trivial and reflect some approximations. The first property is true 
because we follow the dynamics on time scales much larger than the typical correlation time for the 
sources. The isotropy is ‘accidental’ (at least in the relativistic context at hand), in the sense that it 
is specific to the lowest-order approximation at weak coupling and to the ultrarelativistic limit for the 
external particle^ [45]. 

By using a properly discretized version of the stochastic equation (2.14) (see e.g. [45]), one can 
deduce the following evolution equation for the average of the particle momentum squared: 

^ + (2.15) 

For this equation to be consistent with the approach towards the thermal distribution^ fp (x e“^/^ 
(which in turn implies {p) = 3T and (p^) = 12T^), one needs to fulfill the Einstein relation q = 4Tp 
between diffusion and drag. This is indeed guaranteed by the fluctuation-dissipation theorem for 
thermal correlations. 

Let us now assume that the energetic gluon enters the medium at t = 0 with a large momentum 
Po ^ T oriented along the z axis (f = 3) : Vz ^ = 1. At early stages, the longitudinal momentum 
remains large, Pz 3> T, and the effects of fluctuations are unimportant: p± <C Pz and Vz — 1. During 
these stages, one can take the average in Eq. (2.14) to deduce d{pz)/dt ~ —p and therefore {pz{t)) — 
po — pt. This shows that the particle loses most of its energy, from the initial value po T down to 
a value p ^ T where diffusion effects start to be important, over an interval At ~ pg/^ = (Po/^)^rei, 
with 

47^2 

t,ei = ^ ~ a^T\n{l/a) ‘ 

From that moment on, the particle approaches the thermal distribution quite fast, over a time At ~ tj-ei, 
under the combined effect of drag and diffusion. This can be understood from the fact that its 
momentum broadening increases with time like (p^) ~ (3/2)gAt, cf. Eq. (2.15), and hence it becomes 

^In general, {q/2)S'^ in the r.h.s. of the noise correlator gets replaced by = qtv'V -|- {q/2){S'^ — {;*{)■’), where the 
longitudinal (qi) and transverse (q) momentum diffusion coefficients are different from each other. But for a massless 
energetic particle and in the lowest, leading-logarithmic, approximation, it so happens that qi = q/2; see e.g. [45, 47, 49]. 

®The probe gluon is here treated as a classical particle, hence its momentum distribution in thermal equilibrium is 
the relativistic version of the classical Maxwell-Boltzmann distribution, also known as the Maxwell-Jirttner distribution. 
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of 0{T‘^) after a time At ~ /q ~ t^ei- Clearly, when pQ ^ T, the total duration of the thermalization 
process is controlled by the first period — the energy loss via drag — and is of order tth(Po)) with 


(2.17) 

Let us now switch on the branching dynamics, on top of the elastic collisions. From the previous 
subsection, we know that the incoming gluon with energy pQ ^ T has a lifetime At ~ tbr(po) before it 
undergoes a first democratic branching. By comparing Eqs. (2.5) and (2.17), it is clear that tbr(po) ^ 
tth(po) so long as po ^ T. Indeed, 


thvip) ^ J_ / q ^ 1^ 

tth(p) yp’ 


(2.18) 


where we have used q ~ for the weakly-coupled QGP. This implies that the incoming gluon 

disappears via branching before having the time to lose a substantial fraction of its original energy via 
drag. A similar conclusion applies to all the successive generations within the ensuing gluon cascade, so 
long as the respective momenta are hard, p ^ T: the elastic collisions cannot significantly modify the 
kinematics of the hard gluons during the time interval between two successive democratic branchings. 

However, the situation changes in the later stages of the cascade when, as a result of successive 
branchings, the gluons have been degraded to lower energies p ~ T. Then, the various time-scales 
previously introduced become degenerate (at least, parametrically), 

tb,(r) ~ t.h(r) ~ ~ . ( 2 . 19 ) 

meaning that the various processes start to compete with each other. Before they have the time 
to branch again, the gluons with p ~ T can lose a substantial fraction of their energy towards the 
medium and also suffer a considerable broadening of their momentum distribution. Such (drag and 
diffusion) processes will naturally drive the soft gluons towards a thermal distribution in momentum 
which, within the present approximations, appears to be the classical, Maxwell-Boltzmann, distribu¬ 
tion. In reality though, our approximations fail to properly describe the final equilibrium state, for 
the reasons already explained. The gluons from the jet which approach thermal equilibrium cannot 
be distinguished anymore from the gluons in the plasma, so a proper theoretical description of the 
late stages should rather follow the complete gluon distribution. For the latter, the effects of the 
quantum statistics are essential (since the occupation numbers are of order one) and the final distri¬ 
bution in equilibrium must be of the Bose-Einstein type. When this equilibrium distribution is finally 
reached, the branching process naturally stops, because of the compensation between splittings and 
recombination (see the discussion in Sect. 3 below). Thus the cascade effectively ends at a scale ~ T. 

This discussion shows that the characteristic time scale for the thermalization of a mini-jet — the 
duration of the overall process which starts with the emission of a primary gluon with energy po S> T 
and ends up with the thermalization of its soft branching products — is controlled by the branching 
part of the dynamics and hence is of order tbr(po). As discussed around Eq. (2.18), this time scale 
ibr(po) is much smaller than tth(po) — the would-be thermalization time for a gluon with the same 
initial energy in the absence of branchings. This is so since the branchings are more efficient than the 
elastic collisions in degrading the energy of the relatively hard (p S> T) constituents of the jet. 
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In the previous considerations, we have focused on the relaxation of the energy-momentum dis¬ 
tribution, but we ignored the spatial distribution of the gluons inside the jet. In particular, we have 
discarded the fact that this distribution is highly inhomogeneous to start with — the leading particle 
and, more generally, all the energetic constituents of the jet with p ^ T propagate at the speed of light 
and are concentrated near the light-cone (z = t) —, a feature which could prevent, or at least delay, 
thermalization. Also, we have implicitly assumed that the dynamics responsible for thermalization at 
the low-energy end of the cascade {p ~ T) does not affect the branching dynamics at higher energies 
pi$> T. In order to address such complex issues, test our various assumptions, and thus firmly estab¬ 
lish the physical picture that we previously exposed, we need more explicit calculations. A suitable 
formalism in that sense will be described in the next section. 

3 The kinetic equation for the longitudinal dynamics 

In what follows we shall propose a relatively simple kinetic equation which captures the general 
dynamics exposed in the previous section, in the sense that it respects the various parametric estimates 
and thus is expected to reproduce the correct physical picture, and which allows for explicit studies, 
via semi-analytic and numerical techniques. 

The starting point is a kinetic equation introduced in Ref. [36] and thoroughly derived in Refs. [43] 
(see also Refs. [27, 29] for a careful analysis of the quantum branching process, which justifies treating 
the successive branchings as independent from each other), with the schematic structure 

• Va,) /(t,®,p) = Celif] + Cbr[f] ■ (3.1) 

Here, / is the gluon occupation number defined in Eq. (2.7), which a priori refers to gluons from both 
the jet and the surrounding medium®, Cei[/] is a collision integral encoding the effects of 2 —?■ 2 elastic 
collisions, whereas Cbr[/] encodes the inelastic processes, like 2—7-3 collisions, which lead to (collinear) 
branchings. The most general structure of these collision terms can be found in [43, 50], while a simpler 
form of Cbr[/]) which is sufficient for our present purposes, is presented in [36]. These general collision 
terms are non-linear with respect to the gluon occupation number, in such a way to respect the detailed 
balance principle and the quantum statistics. Accordingly, the thermal Bose-Einstein distribution is 
a fixed point for both Cei [/] and Cbr [/] • 

The general equation (3.1) is difficult to solve in practice, because of the non-linear effects alluded 
to above and also because of the generally complicated structure of the two collision integrals, which 
involve multi-dimensional momentum integrations together with complicated kernels (see [43, 47] for 
details). Numerical solutions have been presented for special limits of this equation [39, 51, 52], which 
however do not cover the present physical situation. Here, however, we shall follow a different strategy: 
motivated by the physical discussion in the previous sections, we shall propose a simplified version 
of the kinetic equation, which allows for efficient numerical studies and even for piecewise analytic 
solutions. Similar approximations have been already used in the literature, separately for the elastic 
collisions and for the branching dynamics, but here we shall combine them for the first time in a study 
of the jet thermalization. 

®To Eq. (3.1), one should add corresponding equations for the quark and the antiquark occupation numbers, but here 
these are not necessary, since the (anti)quarks will only appear as constituents of the thermal medium. 
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First, we would like to write an equation for the glnons from the jet alone, with the surronnding 
medium treated as a thermal bath which influences the jet dynamics but is not significantly disturbed 
by the latter. The omission of the back-reaction is indeed well justified, given that the gluonic system 
produced via branching remains dilute (/(p) 1) down to the thermal scale p ~ T, as we have seen. 

By the same token, the kinetic equation can be linearized w.r.t. the occnpation number / for the 
gluons in the jet. 

The distinction between the gluons in the jet and those in the medium is strictly possible for 
the energetic quanta with p ^ T, which control most stages of the branching process, bnt it becomes 
ambiguous at later stages, notably during the thermalization process, where all the gluons have p ^T. 
Notwithstanding, the discnssion in the previous sections shows that the dynamics is rather smoothly 
changing around p ~ T, where the relevant time scales become commensurable with each other, cf. 
Eq. (2.19). Hence, there shonld be no danger (to parametric accuracy, at least) with extrapolating 
down to p ~ T an equation which is strictly valid for p ^ T. With that in mind, one can deduce 
rather simple approximations to the collision integrals in Eq. (3.1). 

Consider the elastic collisions first, which preserve the number of particles. Due to the infrared 
singnlarity of the Coulomb exchanges, the dominant contribution to the elastic collision integral comes 
from the small angle scatterings, i.e. from collisions where the momentum transfer between the 
colliding particles — here, a gluon from the jet and another one from the thermal bath — is much 
smaller than the individual momenta of these particles. More precisely, this property holds to leading 
logarithmic accuracy, since the momentnm exchanges within the range <C 'C T produce a 
Coulomb logarithm ln(T‘^/m‘j^) ~ ln(l/as) in the relevant transport coefficient (see below). Within 
this approximation, Cei[/] can be replaced with the Eokker-Planck dynamics, i.e. the snm of diffusion 
and drag [44-47] 


Ceii/j - \ 



(3.2) 


In writing this expression, we have exploited the isotropy of the diffnsion tensor and the Einstein 
relation between drag and diffusion, as already discnssed in relation with Eq. (2.14). In fact, there is 
a one-to-one correspondence between the Langevin dynamics described by Eq. (2.14) and the Eokker- 
Planck dynamics encoded in the r.h.s. of Eq. (3.2). In particular, it is easy to check that the above 
collision term admits the classical thermal distribution (the Maxwell-Boltzmann distribntion for a 
massless particle) as a fixed point: Cei[/eq] = 0 for feq{p) = where k is independent of p, 

but otherwise arbitrary. Clearly, one cannot expect glnons to obey a classical distribution in thermal 
equilibrium, yet one can rely on this Eokker-Planck dynamics to qualitatively study the approach 
towards equilibrium for the gluons from the jet, which have low occupancy. The relevant mechanism 
at work (elastic collisions with the thermal particles) and the characteristic time scales (namely tth(p) 
for energy loss, cf. Eq. (2.17), and trei for momentum broadening, cf. Eq. (2.16)) are correctly captured 
by Eq. (3.2), to parametric accuracy at least. 

Eor a dilute quark-gluon plasma (QGP) and to the leading logarithmic accuracy of interest, q is 
given by [45, 47, 49, 53] 

g = 87r«2iV,ln(^^) I ^[W/,,(l + /,,)+iV;/,,(l-/,J] 

= asN,Tmlln(^-^^), (3.3) 

\ / 
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where is the leading-order result for the Debye screening mass, that is, 


ml = —asT^ {2N, + Nf) , (3.4) 

Nf is the number of active quark flavors, (fc^ax) ~ is the maximal momentum transfer squared 
between the probe gluon and the medium constituents, and we have included contributions from both 
thermal gluons and thermal quarks, with respective occupation numbers 


/be(p) 


1 


p ^ > 
er — 1 


/fd(p) 


1 

7 ■ 

er -I- 1 


(3.5) 


For the inelastic collision integral Cbr [/] j one can use the corresponding approximation in Ref. [36] , 
which describes medium-induced gluon branching in the LPM regime, with the BDMPSZ branching 
rate shown in Eq. (2.1). The original equation in [36] involves both splitting (1 —)• 2) and recombination 
(2 —7- 1) processes, but only the splitting terms survive after linearizing w.r.t. the gluon occupation 
number. The ensuing expression reads (see also Refs. [27, 29]) 


Chrif] ^ 


thr{p) 


dxlC{x) 


^2 ^ z ^ 


(3.6) 


with tbr(p) as defined in Eq. (2.5) and' 


JC{x) 


[1 — x{l — x)] 2 

3 

[x(l — x)] 2 


(3.7) 


We recall that x and 1 — x represent the longitudinal momentum fractions of the daughter gluons. 
The two terms in the r.h.s. of Eq. (3.6) are recognized as the gain and loss terms associated with a 
collinear splitting: in the gain term, a gluon with 3-momentum p is produced via the splitting of a 
parent gluon with momentum p/x ; in the loss term, a gluon with momentum p disappears because 
it splits. The integrand has singularities at x = 0 and x = 1, which however cancel between the gain 
and loss terms, and the integral is well defined. One can easily check that the branching integral (3.7) 
preserves the total energy: f d^p jpj Cbr[/](p) = 0. 

Notice that the r.h.s. of Eq. (3.6) does not vanish when / approaches the thermal distribution 
feqip) oc (The only fixed point of this particular collision integral is the turbulent spectrum 

/ oc [28, 36, 48].) However, from the previous discussion, we also know that the branching 

term in Eq. (3.6) is strictly correct only when p ^ T. On one hand, the mechanism which triggers 
radiation changes around p ^ T, from multiple scattering to single scattering, so the correct branching 
rate at lower momenta should be of the Bethe-Heitler type. On the other hand, the gluons from the 
jet having p < T cannot be distinguished from the thermal gluons; so, in this soft region of the phase- 
space, the collision integrals should involve the total occupation number (medium plus jet), including 
the associated non-linear effects. This total occupation number efficiently relaxes to the Bose-Einstein 
distribution. When this happens, the non-linear terms — which are omitted in Eq. (3.6), but would be 

^Strictly speaking, the jet quenching parameter q which enters the expression (2.5) for the branching time is not 
exactly the same as that occurring in the Fokker-Planck term, Eq. (3.2), because of the difference between the respective 
transverse scales (fc^ax): one has {k^s.x) ~ for the Fokker-Planck dynamics and, roughly, (fc^ax) ~ 'Jfbr for 

the branching of sufficiently energetic gluons (see the discussion in [54, 55]). Here however we shall ignore this subtle 
difference, which is anyway small, due to the weak, logarithmic, dependence upon (fc^ax), cf. Eq. (3.3). 
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present in a more general version of Cbr valid at soft momenta [36, 43] — will also stop the branching 
process. 

However, this physical mechanism for stopping the branchings is not included in our linear equa¬ 
tion, which applies to the gluon distribution created by the jet alone. To cope with that while keeping 
the formalism as simple as possible, we shall cut off by hand the branching process at some arbitrary 
scale p* ~ T. In practice, we shall enforce the condition p > for all the particles participating in a 
1 —)■ 2 splitting process, that is, we shall require p > p^ for the daughter gluons and hence p > 2p* for 
their parent. This cutoff p* should be viewed as a free parameter of our model: the dependence of our 
predictions upon this scale, which as we shall see is weak so long as p* remains of 0{T), is indicative of 
the error that we have introduced by neglecting the non-linear terms in the (total) occupation number. 

To summarize, our basic kinetic equation reads (with the compact notation fp = f(t,x,p)) 


d \ 1 . 

— -|- • Va; j fp = 


V 


Vp + ^)/, 




tbrip) 


dx fC{x) 


J-f -If' 

s/f 

1X2 ^ ^ 


(3.8) 


where the symbol f denotes the restricted integration over x (see Sect. 5 for details). This equation 

r 

should be solved with the following initial condition at t = 0: 

^3 


f{t = 0 ,x,p) = _^)5(2 )(p^)^ 


(3.9) 


which represents the leading particle propagating along the z axis with energy E ^ T. 

Eq. (3.8) matches our present purposes: it correctly encodes the dynamics of the relatively hard 
constituents of the jet with p » T and, when extrapolated down to p < T, it also describes (at least 
to parametric accuracy) their approach to kinetic equilibrium — that is, the fact that the gluons from 
the jet individually approach a Maxwell-Boltzmann distribution in momentum, via elastic collisions in 
the plasma. On the other hand, the approach to chemical equilibrium — the evolution of the ensemble 
of the complete gluon distribution (jet-|-medium) towards the quantum Bose-Einstein distribution — 
is not encoded in this equation, but merely mimicked in a rather crude way by the lower cutoff p* on 
the branching process. 

Albeit considerably simpler than the original equations, Eq. (3.8) is still too complicated to be 
solved as it stands, including via numerical techniques. A main source of complication is the spatial 
inhomogeneity inherent in our problem, which is very strong to start with, cf. Eq. (3.9), and plays an 
essential role in the subsequent dynamics. In order to keep the salient features of this evolution in a 
numerically tractable way, we shall project Eq. (3.8) along the longitudinal axis and at the same time 
perform approximations based on the separation of scales pz 3> p± between longitudinal and transverse 
momenta. This separation is physically realized so long as p ^ T, which is the regime where Eq. (3.8) 
strictly applies, but is progressively washed out when decreasing the momenta towards T. Yet, this 
approximation correctly keeps trace (to parametric accuracy, once again) of the separation of time 
scales in the problem: indeed, as already explained, the characteristic time scales for branching, 
Eq. (2.5), and for the thermalization of hard particles, Eq. (2.17), are controlled by the longitudinal 
momenta and become degenerate with each other (and with Eq. (2.16)) only when p ^ T. 

Specifically, by integrating Eq. (3.8) over the transverse phase-space, while at the same time 
approximating p 12^ p^ within the definition of Vz, within fbr(p)) and within the condition p > p^:, we 
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finally obtain 


A 

dt 



fi{t,Z,Pz) 


1. d 
4^ dp 2 


d 

dp, T 


fe{t,z,pz) 


+ 


thr{Pz 


dxJC{x) 


fe (t, z, —) - z,pz) 

y/x V X / 2 


(3.10) 


where v, =pz/\pz\ and f^{t,z,pz) is the longitudinal gluon distribution®, 

ft{t,z,pz,) = f d‘^x_id^p±f{t,x,p). (3.11) 

dzdp^ (27r)-^ J 

In Eq. (3.10) it is understood that the partial derivative dp^ = d/dpz commutes with Vz'- dp^{vzfe) = 
Vzdp^fi- (This prescription follows for the limit px ^ Pz- starting with Vz = Pzjp with p = yjp‘1 + 
one obtains dp^Vz = p\/p^ — which is much smaller than the respective natural value 1M-) 

Eq. (3.10) is the equation that we shall explicitly study in what follows, via a combination of 
analytic and numerical methods. To that aim, it is also useful to remind that the branching integral 
above admits the turbulent fixed point fi oc which is expected to control the shape of the 

spectrum at intermediate momenta p^: Pz E. 


4 Semi-analytic studies of the kinetic equation 


As explained in the previous section, the two ‘collision terms’ in the r.h.s. of Eq. (3.10) become 
important in different kinematical regions, which are complementary to each other: the high-energy 
region ai p ^ T for the branching term and, respectively, the low-energy region at p < T for the 
elastic collisions. This distinction makes it possible to separately study their physical consequences 
— at least, at a qualitative level. Namely, one can effectively treat the branching process as a source 
of relatively soft gluons, which get injected into the medium at a scale p^ ^ T and subsequently feel 
the effects of elastic collisions, in the form of drag and diffusion. These considerations motivate the 
following, simplified, version of the kinetic equation (3.10) : 


d d , j.. . 

^ + v-jnt.z,p) 


1 . 

dp 




f{t,z,p) 


+ S{t,z,p;p^). 


(4.1) 


(From now on, we shall omit the subscript i on f as well as the subscript z on longitudinal momenta 
and velocities, to simplify notations.) Notice that, when using this source approximation, we implicitly 
assume that the medium acts as a perfect sink: the branching dynamics at p > p* is not at all affected 
by collisions. For consistency, one must also construct the source S{t,z,p) by assuming an ‘ideal’ 
gluon cascade, with wave turbulence, at p > p*. This source has the general structure 


S{t,z,p;p^:) = 6{t - z)6{p - p^:)T{t,p^:), 


(4.2) 


where r(t,p*) is the flux of gluons at the lower end of the cascade at time t. This flux can be easily 
inferred from the previous studies of the ‘ideal’ branching process [28, 30] and will be presented in 
Sect. 4.2.2 below. Our ultimate purpose in this section is to solve Eq. (4.1) with this particular source. 

^Notice that this quantity ft is not an occupation number by itself (because of the integration over the transverse 
phase-space in Eq. (3.11)), so in practice one can very well have ft>l and still use a linear kinetic equation, provided 
one can justify that the actual occupation number is indeed small. 
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In preparation to that, it will be useful to study a couple of simpler cases, which have a physical 
interest by themselves and for which we will be able to obtain exact solutions in analytic form. We 
shall start by considering in Sect. 4.1 a steady source which propagates at the speed of light. Then, 
in Sect. 4.2.1, we shall construct the exact Green’s function for the differential operator appearing in 
Eq. (4.1) (i.e. for the ultrarelativistic Fokker-Planck equation in 1+1 dimensions). Finally, in Sect. 
4.2.2, we shall use this Green’s function to give a semi-analytic calculation of the gluon distribution 
produced by the physical source. 

4.1 Thermalization for a steady source 

In this subsection we shall present an exact solution for the case where the time-dependence of the 
injection rate r(t,p*) in the r.h.s. of Eq. (4.2) can be neglected. This is a good approximation if one 
is interested in the effects of the collisions over a time interval At which is much smaller than the 
characteristic time scale fbr(F) for the evolution of the source via branchings, but much larger than 
frei (in order for the effects of collisions to be indeed significant); that is, frei < Af <C fbr(-F). During 
this time At, the source can be effectively treated as ‘frozen’ and the corresponding distribution at 
z <t is expected to depend only on t — z. 

For convenience, we choose to normalize the injection rate as r(f,p*) = T. (This brings no loss of 
generality since the equation is linear.) Also, within the context of this subsection, it is preferable to 
denote the energy of the soft gluons as po, rather than p*. With these conventions, Eq. (4.1) becomes 
(below, a prime denotes a derivative w.r.t. p) 

+ / = if + vfY + 5{i - z)6{p - Po), (4.3) 

in terms of dimensionless variables which measure the respective quantities in natural units, that 
is, in units of frei for all the time and length scales (t = f/freP z = z/tj-ei) and in units of T for 
the various momenta [p = p/T, etc). In what follows, we shall use such reduced variables in most 
formulm, to simplify the notation, but we shall restore the physical units when discussing the physical 
interpretation of the results. Also, we shall drop the hat on the reduced variables (e.g. p —)• p), as the 
distinction should be clear from the context. 

We search for a stationary distribution /(x“,p;po) with x~ = t — z. This function obeys 
^ = if + f)'+ 6ix~)5{p - Po) for p > 0, 

< o (4.4) 

= (/'-/)' forp<0 . 

together with the condition for particle number conservation at p = 0: 

(/' + /)+„* = {/'-/)+„- =<5(x-). (4.5) 

For p > 0 the solution to Eq. (4.4) is found in the form 

fix~,p;po) =/j(p,po)(5(x") + C+(a:")e“P (4.6) 

where /j(p,Po) is the ‘jet front function’ (see also Fig. 2) 

hiP,Po)^ e~P ieP° - l)9ip-Po) + il - e~P) eipo-p), (4.7) 
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Figure 2. The jet front function /j(p,_po) is displayed as a function of p for po/T = 1,5 and 15. 

and C~^{x~) is an unknown function, to be later determined. For p < 0, it is convenient to use the 
Laplace transform of the solution, fs{p) = dx“e“*^ f{x~,P',Po)- By taking the Laplace transform 

in the second equation (4.4), one obtains 

/,(p) = C-(s)e2(^/l^+l) 


for p < 0 . 


(4.8) 


By imposing the conservation condition (4.5) (which actually introduce two constraints), we can 
determine both ‘coefficient’ functions C~{s) and C''’'(s) (the Laplace transform of C~^{x~)) : 

2 


C+(5) = C7-(s) = 


\/8^T1- 1 

After also performing the inverse Laplace transformation, we hnally obtain 


(4.9) 


f{t-z,p;po) =< 


fj{p,Po)S{t -z) + 


ierf 


2 V 2 


+ 1 -I- _B_ 

4 


4^ 


erf 


2p+t-z \ I 1 
2y/2y/t^) 


+ 


i-2p+t-zy 

8(t-z) 

y/^y/t — Z 


e P for p > 0, 


for p <0. 


(4.10) 


We have introduced here the error function 
9 

eif{x) = / dte“*^. (4-11) 

JO 

By inspection of the solution in Eq. (4.10), one can recognize a front which propagates at the 
speed of light with the prohle in p shown in Fig. 2, and a tail at z < t which is localized around p = 0. 
The overall distribution is illustrated in Fig. 3 for the case po = T, which is the most interesting 
one for the physical problem at hand (recall that po = p* is the infrared end of the gluon cascade). 
However, in order to better appreciate the physical content of this stationary distribution, it is useful 
to consider first the high-energy case po T. Then, as visible in Fig. 2, the front profile becomes a 
0-function with support at 0 < p < po. This can be understood as follows: a particle injected by the 
source at time to with po S> T loses energy towards the medium at a constant rate, via drag (recall 
the discussion following Eq. (2.14)); hence, its energy decreases with time according to 


p{t) =P0- (r/trel)(t - to). 


(4.12) 
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Figure 3. The distribution (4.10) as produced by a steady source is shown as a function of t — z and p for 
Po = T. This figure exhibits a front moving along the light-cone at z = t and a thermalized tail at z t — t^ei- 
In displaying the front, the d-function in Eq. (4.10) is regulated as 6e{t — z) = ^ with e = 0.1. 

So long as this energy remains much larger than T, which is indeed the case during a large interval 
t—to — (po/7")irei 3> treh the diffusion effects are negligible and the distribution created by this particle 
can be as well studied by neglecting the second-order derivative / in Eq. (4.3). The corresponding 
solution is easily found as 

/(t — z,p;po) = 6{t — z)9{pq — p)9{p) (without diffusion), (4-13) 

which is indeed very similar to the ‘front’ piece of Eq. (4.10) in the case pq ;§> T. Hence, the ‘front’ 
is built with those particles that have been recently injected by the source, within a time interval 
= {po/T)tre\ prior to the time t of measurement, and which have a still a relatively large energy 
p > T at time t. All the other particles, that have been injected at earlier times t' < t — {po/T)t^ei, 
have been degraded by the viscous drag to energies p '^T, where the diffusion effects are important. 
This becomes clear by inspection of the function fj{p,po) for pQ = 1 \n Eig. 2. 

As a consequence of the competition between diffusion and drag, the gluons with \p\ < T can have 
both positive and negative velocities, hence their distribution moves at a slower speed |u| < 1 This 
explains the depletion visible in /j at p < T (for any po) and also the formation of the tail. At points 
sufficiently far away from the front, such that t — zf$> treh the distribution reaches thermal equilibrium, 
since this is the fixed point of the Fokker-Plank dynamics. Indeed, the solution in Eq. (4.10) implies 

f{t-z,p;po) ^e"IPhr when t - z > t^ei • (4.14) 

To understand the energy balance between the jet and the medium, notice that the external source in 
Eq. (4.3) inserts energy at a rate dEg/dt = poT, whereas the energy carried by the thermalized tail 
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t — z ^ tj-ei increases at a rate d-Ether/dt = T^- (The total energy in the front in independent of time, 
Ej = J dppfj{p,po) ~ treiPoi ™d represents only a negligible fraction of the total energy injected by 
the source over large times.) If po ^ T, then the insertion rate is much larger than the thermalization 
rate, meaning most of the energy is lost via viscous drag — meaning it is transferred via collisions to 
the medium constituents. By choosing po ~ T, i.e. by inserting the particles directly at the medium 
scale, we can minimize the effects of the drag and thus recover most of the injected energy in the 
thermalized tail. 

The qualitative features that we have discovered in this simple example are in fact generic and 
will be recovered in the more general situations to be studied later on. In particular, the peculiar 
structure of the distribution visible in Fig. 3, with a jet localized on the light-cone {z = t) and a 
thermalized tail well behind it {z < t — trei); will also show up for a physical jet initiated by a leading 
particle with energy E ^ T, at least for sufficiently small times t < tbr(-®)- 


4.2 Jet quenching in the source approximation 

In this subsection we shall use the source approximation, cf. Eq. (4.1), in order to unveil generic 
features of the jet evolution in the presence of both branchings and elastic collisions. The solution to 
Eq. (4.1) corresponding to a general source S{t, z,p) can be written as 


f{t,z,p) = 


j dpodzo fcit, z - zo,p,po)fQ{zo,po) 


+ / dpodzo / dt'fGit-t',z-zo,p]Po)S{t',zo,po), 


(4.15) 


where we have chosen the initial condition f{0,z,p) = fo{z,p) and fc is the appropriate Green’s 
function, that is, the solution to the homogeneous equation 


d d \ J. / N 


A 

dp 



fG{t,Z,p) , 


(4.16) 


with initial condition fdO, z,p]po) = 5{z)5{p — po) with pQ > 0. An analytic form for this Green’s 
function will be constructed in the next subsection and then applied to the source representing an 
ideal branching process, in Sect. 4.2.2. 


4.2.1 The Fokker-Planck Green’s function 

The Green’s function for the longitudinal Eokker-Planck equation can be constructed via a mathe¬ 
matical method similar to that described in the previous subsection for the case of a steady source. 
In what follows, we shall omit the details but merely show the starting point equations, which replace 
the previous equations (4.4) and (4.5). After performing Laplace and Eourier transforms with respect 
to t and z respectively, we deduce 

(s -F iQ) - 5{p - Po) = fjg + f'^Q for p > 0, 

< (4.17) 

_ sfsQ - iQfsQ = fsQ - fsQ for p<0. 

together with the following condition for the number conservation 

{f'sQ + fsQ)\p^Q+ = {fsQ - fsQ)\p^Q- , (4.18) 
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where fsQ is a compact notation for a function of two arguments, s and Q, defined as 

POO P 

fsQ^J dte~^^ J (4.19) 

After lengthy but straightforward mathematical manipulations, one finally obtains 


fG{t,z,p;po) = 0{p)fQ{t,z,p;po) + 9{-p)fQ{t,z,p-,po), 


(4.20) 


where (with t > \z\) 


/G(i,^,P;Po) = 


PQ-P t 
e 2 4 

2\/7rt 


(p-Po) 


(p+Po) 


e 4t — e 4t 


d(t — z) 


(p+PO-z) 

4t P 


2^ f P + PO 


t V ^ + -2; 


- 1 


and 


+ 8 V7rt5/2 + 2)-ip + po- zf] erfc 

I it + z){p + po + t-z) 


^ , , (t + z)(p(} + t — z) — p(t — z) _ PZ _ : P+P0 _20_i 

fQ{t,Z,p-,Po) = - - . - -e 2(t-z)+ 2 2(t+z) 4 


47rt^\/t^ — z'^ 

C* r 01 

+ 8V^t5/2 + 2) - (p + po - ^) ] erfc 


(P+PQ-Z) 


-A/t-- 


/ J?o 




t + t — z 


- 1 


(4.21) 


(4.22) 


In these formulas, we have introduced the complementary error function erfc(a:) = 1 — erf(x) (cf. 
Eq. (4.11)). Also, we have used reduced variables p —)■ p/T, t —)• t/trei etc., to simplify writing (recall 
the discussion after Eq. (4.3)). The above expression for fc is normalized to unity w.r.t. these reduced 
variables: f dzdp fdt, z,p;po} = 1. In order to obtain the properly normalized Green’s function in 
physical units, one must divide the result in Eqs. (4.21)-(4.22) by the dimensionless product Ttrei and 
replace all the reduced variables by their physical counterparts (p —)■ p/T, etc). 

The above expression for fc looks quite involved. In order to unveil its physical content, it is 
useful to first consider its simpler version obtained after integrating over z : 


fGit,p;po) = < 


4® 


-p 


erfc 


- e^ erfc 
4 


P + Po-t 

2yft 

-p + Po-t 

2Vi 


+ 


1 


2\/Tri 
e 


(p-pn+*)'^ 
e 4t 


(-p + PQ-f)^ 

e 4t 


2\/Tri 


for p > 0 
for p < 0. 


(4.23) 


This function describes the relaxation of an initial perturbation which is homogeneous in z but localized 
in momentum: /g(0,p;Po) = S{p — po)- Eor po ^ T and sufficiently small times, such that ^ ^5 
/G(t,p;po) is dominated by its Gaussian component at p > 0, that is. 


/G(t,p;po) ^ 


1 (p-PO+t)^ 

—= e 4t 
2\/TTi 


when t <C Po • 


(4.24) 


This describes the damping of the original energy via drag and also the broadening of the longitudinal 
momentum distribution due to diffusion, in agreement with the discussion in Sect. 2.2: 


(p(t)) ~ Po - iT/t,e\)t , (p^) - {pf ~ 2qt. 


(4.25) 
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Figure 4. The ‘reduced’ Green’s function in Eq. (4.23) plotted as a function of p for po = 5 and 6 successive 
values of time: t = 0.1, 1, 2, 3, 5, 10. 


But already for such small values of time, there is a second component (represented by the two 
terms proportional to the complementary error function) which starts growing around p = 0. This 
corresponds to particles which have essentially lost their original energy and pile up around p = 0, 
via diffusion. For larger times t > {po/T)t^^\, this component becomes the dominant one and rapidly 
approaches the thermal distribution (notice that erfc(x) —?■ 2 when x —?• —oo) : 


fG{t,p;Po) - IpI 


when t Po . 


(4.26) 


Turning now to the general case with z-dependence, Eqs. (4.21)-(4.22) show that, at early times, 
t <C {po/T)t^eh there is a remnant of the original perturbation — the ‘jet front’ localized at z = t, 
as described by the first term in the r.h.s. of Eq. (4.21) — whose momentum distribution is however 
degrading with time, in the same way as in Eq. (4.24). With increasing time, this ‘jet front’ is gradually 
washed out and a new distribution develops around p = 0, which is slowly varying in z (for \z\ <C t at 
least). For sufficiently large times, t i$> (po/T)trei, and sufficiently far behind the front, z < t — t^ei, 
this new distribution is thermal in p : 


e-bl 

fGit,Z,p-,po) ~ 


(z-pn)^ 

e 


21/71^ 


when t Po ^ 1 and |z| <C t. 


(4.27) 


This is recognized as the product between the Maxwell-Boltzmann distribution in momentum, cf. 
Eq. (4.26), and the one-dimensional heat kernel describing diffusion in z. Eq. (4.27) shows that 
the large-time distribution is centered around z = {po/T)t^f.\ (the maximal distance travelled by the 
original perturbation until it has lost all its energy, due to drag) and that its longitudinal extent 
grows with time, according to Az(t) ~ Thus, remarkably, the momentum-space diffusion 

originally encoded in the Eokker-Planck equation has also generated spatial diffusion. One may chose 
as a criterion for having a quasi-homogeneous distribution the condition that the longitudinal extent 
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Figure 5. The Green’s function in Eqs. (4.21)-(4.22) plotted as a function ofp and z for po = 5 and four values 
of t: upper line, left: t = 0.1; upper line, right: t = 1; lower line, left: t = 5; lower line, right: t = 20. Note that 
the vertical scales and also the ranges in p and z can significantly differ from one figure to another. For t > 1, 
one can see the emergence and growth of a thermal tail at \z\ t, which for t ^ po is the Gaussian shown in 
Eq. (4.27). 


Az{t) be larger than the location (po/r)trei of the center. This happens for 

(4.28) 

a time scale considerably larger than that required by the thermalization of the momentum distri¬ 
bution. Note finally that by integrating Eq. (4.27) over z and over p we recover, to the accuracy of 
interest, the normalization of the initial perturbation, i.e. f dzdp felt, z,p} = 1, which confirms that 
the whole perturbation has thermalized. If on the other hand one computes the total energy contained 
in the thermalized distribution (4.27) one finds, clearly, f dzdp\p\ fG{t-,z,p) ~ T for t » (po/r)irei- 
That is, out of the total energy po > ^ of the original perturbation, a fraction T/po is carried at large 
times by the thermalized distribution, whereas the remaining fraction (po ~ '^)/Po has been transmit¬ 
ted to the medium, via drag. This conclusion on the energy loss is similar to our previous findings for 
the case of a stationary source in Sect. 4.1. 
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4.2.2 A physical source generated by the branching process 

In this subsection we apply the Green’s function method to the main physical problem of interest, 
namely a source generated by a branching process. More precisely, we consider an ideal branching 
process, for which the splitting dynamics at p > p* is not at all influenced by elastic collisions: the 
medium solely acts as a ‘perfect sink’ which absorbs the energy of the gluon cascade at the ‘infrared’ 
scale p* ~ T. (A more general situation will be studied in Sect. 5.) Under these circumstances, the 
source has the structure shown in Eq. (4.2) with r(t,p*) = J^(E,p*,t)/p*. Here, is the 

energy flux at p* generated at time t by a cascade that was initiated at t = 0 by a leading particle 
with initial energy E ^ T (recall the discussion in Sect. 2.1). This source truly represents a bunch of 
relatively soft particles, which carry all the same energy p*, move together at the speed of light, and 
whose number is evolving in time due to the branching dynamics. 

Within the context of the ideal cascade, one was able to obtain an exact analytic result for this 
function E{E,p^,t) [28]. Strictly speaking, the analysis in [28] required two additional assumptions, 
which are not essential from the viewpoint of physics, but simplify the mathematical manipulations: 

(a) the energy E of the LP is not too large, namely it obeys® E < uJc{L) = where the upper 
limit 0 Jc{L) = Wbr(U)/Q;^ is parametrically larger than a;br(T) at weak coupling; 

(b) the kernel (3.7) within the BDMPSZ splitting rate is replaced by its simplified version JCq{x) = 

3 

l/[x(l — x )]2 , which preserves the correct behavior near the singular endpoints at x = 0 and x = 1. 

Under these assumptions, the energy flux is obtained as [28] 


E{E,p^,t) 


2ttE 


tliE) 


dt 


flow ' 


(4.29) 


where tbr(U) is the branching time introduced in Eq. (2.5) (the typical time after which a parton with 
energy E undergoes the first democratic branching) and 

AE{iowiE,p^,t) = E[l- = E[l- 


with t<;br(^) = di^qt^i is the energy which accumulates into the soft modes with p < P* after a time t. 
The above results strictly apply for p* <C Ci;br(^)) or, equivalently t ^ tbr(p*) ~ tveh arid within that 
regime they are independent of p*. In fact, in the absence of the sink at p*, this whole energy would 
accumulate in a condensate at p = 0 [28] . 

Eq. (4.30) confirms that the time scale tbr(711) plays the role of the lifetime of the leading particle 
w.r.t. democratic branchings. At small times t <C t\,j;{E), one can expand the exponential there to 
lowest order and thus find 


AEfiowit) — 7rwbr(i) when t <C t]:,^{E) . (4.31) 

Alternatively, this estimate is correct for a given time t provided the energy E of the LP is sufficiently 
high, E S> ujhr{t)- This result is in agreement with (2.6) up the replacement v — )• 27r, due to the use of 
the approximate kernel /Co(x). It shows that, at small times, the energy loss via flow is independent 
of E (essentially, because the LP does not ‘feel’ the change in its own energy due to flow) and that it 
grows with time like . As explained in Sect. 2, this ‘small time’ (or ‘high energy’) regime is the most 
interesting one for jets at the LHC, where one indeed has E ^ cubr(L), with L the size of the medium. 

®The generalization of the subsequent results to more energetic jets with E ^ uJc{L) can be found in Ref. [30]. 
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Figure 6. The distribution (4.32) produced by the source in Eq. (4.29) with thr{E) = 6 and p* = 1 is plotted as 
a function of p and z for four values of t: upper line, left: t = 2; upper line, right: t = 5; lower line, left: t = 8; 
lower line, right: t = 20. For relatively small times t < ti^,.{E), the source is still active and the distribution 
is thermal only in the tail at z < t — trei- For larger times t > thr{E) +p*, the source has essentially decayed 
and the distribution is thermal at any z. For t = 20 there is no significant difference between the exact result 
displayed above and the respective prediction of the diffusion approximation (4.33). 


For larger times, such that t ~ tbr(-E') or, equivalently, E ~ cubr(i), the LP disappears via demo¬ 
cratic branching and its whole initial energy is carried away by the flow. This is indeed consistent 
with Eq. (4.31) which shows that AEflow(t) — E when t > tbr(-E')- From a physical viewpoint, this 
‘large time’ (or ‘low energy’) regime better corresponds to the primary gluons radiated by the LP, 
which evolve into ‘mini-jets’. 

It is also interesting to notice the time dependence of the energy flux in Eq. (4.29): this rises 
linearly with t at small times, then reaches a maximal value around t = tbr(-E')) and rapidly vanishes 
at larger times. This means that the production rate for soft gluons is largest towards the late stages 
of the branching process, i.e. for t ~ tbj-(E), just before the LP dies away. 

We now return to the solution to the Fokker-Planck equation for the physical source at hand. 
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Using Eq. (4.2), the integrals over zq and p* in Eq. (4.15) can be immediately performed to yield 


f{t,z,p) = — [ dt'fG{t-t',z-t',p-,p^)T{t'). (4.32) 

P* Jo 

It seems difficult to analytically perform the remaining integral over t', but this can be numerically 
computed, with the results shown in Fig. 6 (in terms of reduced variables^*^). These results can be 
understood as follows: For relatively small times, t <C tbr(-£')j the source (the leading particle) is 
still present and the gluon distribution is quite similar to that produced by a stationary source, as 
shown in Fig. 3: it exhibits a front at z = t and 0 < p < p^ (but with the edges smeared out by 
diffusion), which represents the gluons that have been recently emitted, and with a tail at z < t 
and peaked in momentum at p = 0, which describes the gluons which have experienced the effects 
of collisions. For later times t > toriE), the source has disappeared via democratic branching and 
the gluon distribution looks quite similar to that produced by a localized source at late times, cf. 
Fig. 5: a thermal distribution in momentum which extends in 2 via diffusion and which peaks around 
2 = t\yr{E) (the maximal displacement of the source before it dies away). 

In fact, for sufficiently large time, t — t^riE) ^ {p^/T)tT-ei ~ t^ei, one can use the diffusion 
approximation for the Green’s function, Eq. (4.27), to deduce (in reduced variables, cf. footnote 10) 


^,2 


f{t, z,p) ~ dvra 


:.-bl 


dt't' 


■ 4{t-t') 




P* Jo y/d7r{t — t') 


(4.33) 


This represents the medium perturbation that would be left over by a relatively soft mini-jet {E <C 
after it thermalizes. By integrating this late-time distribution over z and p, it is easy to check 
that it encompasses all the particles generated by the source, f dzdpf(t,z,p) = E/p^,, as it should, 
and that is contains a fraction T/p^ of the total energy: 


■^-^ther 


j dzdplpl f{t,z,p) 


large time 





when t > tbr(-E) • 


(4.34) 


We thus conclude that by choosing p^, =T one can ensure that the whole initial energy of the source is 
eventually recovered in the thermalized gluon distribution: the energy loss via viscous drag is negligible 
since the gluons are directly injected at the thermal scale. 


5 Numerical studies of the kinetic equation 

The general discussion and the parametric estimates for the energy loss presented in Sect. 2, as 
well as the explicit calculations using the Green’s function method in Sect. 4.2.2, were based on an 
important physical assumption: the fact that the gluon cascade generated via multiple branchings is 
not modified by the elastic collisions responsible for thermalization and hence it can be modeled as an 
ideal branching process, along the lines of Refs. [28, 30] — that is, a turbulent cascade, for which the 
medium acts as a perfect sink at the lower end {p ~ T) of the cascade. The validity of this assumption 
is far from being obvious, as shown by the following argument: via elastic collisions, the soft gluons 

When rewriting Eq. (4.1) in terms of reduced variables, we find that the dimensionless version of the source in 
Eq. (4.2) reads (we restore the hat on reduced quantities, for more clarity): S = 6{t)S{t — z)5{p — Pf){J-/p*), where J- = 
TjT'^ = STva^iexpl—TvP/t^j.{E)} depends upon the energy E only via the reduced branching time ihr{E) = tbr(E)/trei. 
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are redistributed in phase space, in such way to match a thermal distribution, and for p ^ T the latter 
is qnite different from the scaling spectrum produced by the turbulent cascade. 

In this section, we shall give up the ‘perfect sink’ assumption and present a detailed nnmerical 
study based on the kinetic equation (3.10). This equation too is a rather simplified version of the 
actnal dynamics, as explained in Sect. 3, but as compared to the source approximation in Sect. 4 
it has the merit to include an explicit infrared cutoff p* ~ T in the branching process and also the 
interplay between branching and thermalization at p > p*. 

5.1 Setting-up the problem 

As in the previous section, it is convenient in practice to measure all the momenta in units of T and 
all the space-time scales in units of trei, that is, to nse the reduced variables introduced in Eq. (4.3). 
In terms of these variables, the kinetic equation (3.10) reads (with v = p/|p|) 

{dt + vdz) f{t,z,p) = dp{dp + v)f{t, z,p) 

H- —Y [ dxlC{x) 

tbr(T)p^ J 

where we recall that the subscript r on the integral over x indicates the condition that both daughter 
gluons in a splitting process be harder than the ‘infrared’ scale p* ~ T (the lower end of the cascade). 
The precise kinematical conditions are as follows: for the gain term, p > p* and p(l — x)/x > p^, 
whereas for the loss term, xp > p* and (1 — x)p > p*. For definiteness, we chose this scale p* to be 
exactly equal to T (i.e. p* = 1 in Eq. (5.1)). 

The ratio trei/^br(^) which appears in front of the branching term in Eq. (5.1) is parametrically 
of order one for the weakly-coupled quark-gluon plasma. For what follows, it is convenient to choose 
this ratio to be exactly eqnal to one. This choice is also reasonable from a physics standpoint: it 
corresponds to the condition q = (or, equivalently, = 1), which is satisfied by the 

following values for the physical parameters: 

a = 0.3, T = 0.5 GeV, ^ = 1 GeV^/fm ~ 0.2 GeV^ , Gel = 1 fm, (5.2) 

which are indeed consistent with the current phenomenology. 

The equation thus obtained will be solved numerically, with the initial condition 

fit = 0, 2 , p) = 5(p - E)5iz) ^ ^ ( 5 . 3 ) 

TT 

where the product of d-functions is regulated as shown in the r.h.s. The picture that we expect in the 
light of the general discnssion in Sect. 2 is as follows: 

(i) For sufficiently small times t <C Gr(.£■)) the leading particle should survive and carry most of 
the total energy. The energy lost towards the medium should be comparatively small and follow the 
law shown in Eq. (2.6) ; that is, it should be of order Wbr(t) and thus grow with t as 

(ii) For larger times t > tbr (£')5 the LP should disappear via democratic branching and the 
energy loss shonld be of the order of the total energy E. 

For what follows, one should keep in mind that some of the assumptions nnderlying this picture 
might not be well satisfied when solving Eq. (5.1) in practice. For instance, in Sect. 2 we have assumed 




(5.1) 
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the medium to act as a perfect sink for the energy carried away by the branching process, which in turn 
allowed for a well-developed phenomenon of turbulence. This is an important hypothesis, implicitly 
assumed in the previous literature, for which our subsequent study will provide an explicit test. 

5.2 The gluon spectrum 

Before we describe the full picture of the gluon cascade in the longitudinal phase-space {z,p), let us 
present the results for the distribution integrated over z, that is, the gluon spectrum 

^ = j dzf{t,z,p). (5.4) 

Clearly, this function can be obtained by solving directly the homogeneous (in the sense of independent 
of z) version of Eq. (5.1), with initial condition f{0,p) = 6{p — E). This is a relatively simple numerical 
problem, which in particular allows for an extensive study of the role of the lower cutoff at p = p* on 
the branching process. To that aim, we shall consider three cases: (a) an ‘ideal’ branching process, 
which involves no infrared cutoff and develops a clear phenomenon of turbulence (the corresponding 
equation is obtained by keeping only the branching term in the r.h.s. of Eq. (5.1) and letting p* = 0); 
(b) a branching process with a sharp infrared cutoff at p* <C E, as described by Eq. (5.1) without the 
Eokker-Planck terms, and (c) the complete dynamics (branchings with infrared cutoff p* and elastic 
collisions), as described by the homogeneous version of Eq. (5.1). 

As already mentioned, the ideal branching process has been extensively studied in the literature 
and, in particular, exact analytic solutions have been obtained [28] for the case of a simplified kernel 

3 

IC{x) —)• JCo{x) = l/[x(l — x)] 2 . In our numerical study, we use the full kernel in Eq. (3.7), but the 
solution is qualitatively similar to that presented in Ref. [28]. Namely, for sufficiently small momenta 
p <C E, the spectrum exhibits the scaling law f{t,p) oc 1/p^/^, which is a fixed point of the branching 
kernel and the signature of wave turbulence. Eor not too large times t <C tbr(E), the leading particle is 
visible in the spectrum, as a pronounced peak just below p = E. The width of this peak increases with 
time (due to radiation) and eventual becomes of order one — meaning that the LP undergoes its first 
democratic branching — when t ~ tbr(E). Eor even larger times, the scaling law ~ 1/p^/^ is still visible 
at small p, but the spectrum is suppressed as a whole, since the energy flows via multiple branching 
and accumulates at p = 0. When t ^ tbr(E), the whole energy E ends up in this ‘condensate’. This 
behavior is clearly visible in the numerical results displayed in Eig. 7. 

After introducing the infrared cutoff p*, the gluons with p < 2p* cannot split anymore, as there 
is no phase-space available to the daughter gluons. Hence, instead of falling at p = 0, gluons start 
accumulating in the bins at p > p*. As a result, the spectrum above p* deviates from the scaling 
spectrum; it shows an excess (‘pile-up’), which is particularly marked at p* < p < 2p*, where it looks 
like a bump. Both the size of this excess and its extent in p above p* are increasing with time, as 
clearly visible in Eig. 7. This can be understood as follows: a gluon with, say, p = 3p* has more 
chances to be created via the decay of parent gluons with p^ p^, (for which the kinematical constraint 
is relatively unimportant) than to disappear via a decay (since a significant fraction of the phase-space 
for its decay, that at p < p*, is not accessible anymore). 

On physical grounds, it is quite clear that this pile-up cannot be entirely physical: gluons with 
p ~ r can efficiently lose energy towards the medium via elastic collisions and hence they should fall 
into the bins at lower energies |p| < T. We thus expect the pile-up to be considerably reduced and 
possibly washed out after also including the elastic collisions, as represented by the Eokker-Planck 
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t = 0.02 V(E) = 0.2 t,e| 



t = 0.1 tbr{E) = 0.9 t,e| 



t = 0.5 tt,r(E) = 4.7 Vei t = ^(E) = 9.5 




Figure 7. The gluon spectrum f{t,p) for a spatially homogeneous distribution, for an initial energy E = 90T 
and 4 values of time. The three curves correspond to (a) an ideal branching process (the red, continuous, 
line), (b) a branching process with infrared cutoff p* = T (orange, dotted line), and (c) the full process with 
branchings and elastic collisions (green, dashed-dotted line). On the vertical axis, the spectrum is multiplied 
by (p/T)^/^ to render manifest the scaling behavior for the ideal branching process. 


terms in the r.h.s. of Eq. (5.1). This expectation is conhrmed by the numerical solution to the 
homogeneous version of Eq. (5.1) (whose results are shown too in Eig. 7), but only partially: the pile- 
up in the spectrum is indeed reduced by the elastic collisions, but the deviation with respect to the 
scaling spectrum remains quite large — so large, that there seems to be no scaling window in practice. 
If true, the last conclusion would also imply that the physical results are very sensitive to the details 
of the mechanism which stops the branching process and which in our analysis has been only crudely 
mimicked by the infrared cutoff p*. Fortunately though, these last conclusions are not fully right and 
the numerical results exhibited in Fig. 7 are in this respect quite misleading. The gluons which appear 
to accumulate on top of the scaling spectrum in this figure are actually located at different values of 
z. These are relatively soft gluons, which undergo strong diffusion as a consequence of collisions and 
thus separate from each other and also from the more energetic constituents of the jet (which keep 
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propagating along the light-cone at z = t). Hence the ‘pile-up’ visible in the curves denoted as ‘full’ in 
Fig. 7 is merely an artifact of integrating the gluon distribution over the longitudinal coordinate z: it 
comes from the superposition of a nearly ideal branching spectrum in the front of the jet at z ~ t and 
of nearly thermal spectra in the tail of the jet at z <C t. This will be demonstrated by the subsequent 
analysis, where the z-distribution is kept explicit. 

5.3 Jet evolution in longitudinal phase-space 

In this subsection we present numerical solutions to the complete equation (5.1) with initial conditions 
of the type shown in Eq. (5.3). We consider two values for the initial energy, E = 25T and E = 90T, 
which correspond to rather distinct physical situations. The first value E = 25 T (= 12.5 GeV 
according to Eq. (5.2)) is quite low and is representative for a mini-jet radiated by a leading particle 
with a much higher energy Eq > 100 GeV. This is one of the typical mini-jets which control the energy 
loss in the case where the LP crosses the medium along a distance L ~ tbr(-E') = 5trei = 5 fm. The 
second value E = 90 T = 45 GeV is closer to the energy of an actual jet at the LHC and in particular 
is large enough to ensure that the respective leading particle does not disappear into the medium; 
indeed, the respective branching time tQ^{E) ~ 9.5 fm is larger than the typical distance L < 8 fm 
that the LP might travel across the medium in the experimental situation at the LHC. 

5.3.1 The gluon distribution and the energy density 

The general features of the evolution of the gluon cascade produced by a high-energy jet can be 
appreciated by inspection of Figs. 8 and 9, which show the phase-space distribution of the gluon 
number f{t,z,p) and, respectively, the gluon energy \p\f{t,z,p), for the two energies of the LP, 
E = 90T and 25 T, and four values of time: t/trei = 0.95, 4.7, 9.5, and 14. These particular values for 
t have been chosen since, with our present conventions (i.e. tbr(-®)/^rei = \/E/r), they correspond to 
rather special values for the case of a jet with E = 90T] namely, they amount to t/tbr(90) ~ 0.1, 0.5, 1, 
and 1.5, respectively. Since the natural time scale for the jet evolution is tbr(-F), let us also list the 
corresponding values for the softer jet with E = 25T : one roughly has t/tbr(25) ~ 0.2, 1, 2, and 3, 
respectively. To gain more intuition about these time scales in physical units, it is useful to recall that 
trel = 1 fm for the medium parameters in Eq. (5.2). 

When discussing Figs. 8 and 9, it is natural to group together those plots which correspond to 
different values of the energy E, but similar values of t/tQ^{E), because they refer to similar stages in 
the evolution of the jet via branching. But even for identical values of t/t\,j-{E), one should still expect 
some differences between the two cases, E = 90 T and E = 25 T, because the physics of thermalization 
introduces an additional energy scale in the problem — the infrared cutoff p*. 

Consider first Figs. 8 (a)-(c) and Figs. 9 (a)-(c), which illustrate the evolution at early stages, 
t < tbr(-F)- Figs. 8 (a)-(c) show that, already for such early times, most of the particles are relatively 
soft (p ~ p* = T), meaning that they are products of radiation. In particular, those particles which 
at a given time t have an energy |p| smaller than a;br(i) = cx^qt"^, are generally produced via multiple 
branchings, that is, they belong to gluon cascades generated via democratic branchings by primary 
gluons with p ~ Wbr(i)- But so long as t <C thr{E), most of the energy is still carried by the LP, as 
manifest in Figs. 9 (a)-(b): the energy distribution is peaked at a value which is smaller than, but 
comparable to, the original energy E. When t^ei ^ t tQj;{E), the energy loss by the LP and its 
longitudinal broadening are both controlled by soft branchings and hence they grow with time like t^. 
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Figure 8. The time evolution of the phase-space distribution f{t,z,p) produced by a jet with initial energy 
E = 90 T (left) and respectively E = 2bT (right), plotted for exactly the same values of time. 
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Figure 9. The time evolution of the phase-space energy density \p\f/T, for the same conditions as in Fig. 8. 
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Figure 10. The gluon distribution f{t, z,p) for E = 25T is represented as a function of p at time t = 2t^e\ and 
for 3 different values of z : z = 0, z = Eei, and z = 2<rei- For z < t — Eeu the momentum distribution is nearly 
thermal, but with an overall strength which depends upon z : f{z,p) ~ /o(z)e“I^F'r_ ^ = t on the other 

hand, the distribution is far from thermal equilibrium, albeit it is still strongly peaked near p = T. Similar 
conclusions hold for the distribution at later times and for a jet with E = 90T. 


Fig. 9 (c) shows another interesting feature: for t ~ 0.5 one sees a second peak emerging 
in the energy distribution at p ~ T. This demonstrates the strong accumulation of gluons towards 
the lower end of the spectrum which in turn reflects a limitation of the medium capacity to act as a 
‘perfect sink’. We shall later return to a more detailed study of this phenomenon (see notably Figs. 11 
and 14 and the associated discussions). 

As also visible in Figs. 8 (a)-(c), the approach to thermalization in the tail of the distribution 
at z < t is noticeable already at such early times t < thr{E). There is indeed a substantial number 
of gluons which remain behind the LP (i.e., which do not travel at the speed of light) and whose 
momentum distribution is nearly thermal. This is more clearly illustrated by the plots in Fig. 10, 
corresponding to E = 25 T, which show the momentum distribution at t = 2 Eel = 0.4fbr(25) and for 
different values of z: the shape of this distribution is close to the exponential at any z < t — Eei- 

Consider now later times t > E^{E), where one expects the LP to disappear via democratic 
branching. Figs. 9 (d) and (e) confirm that, when t ~ Ej.{E), there is no visible trace of the LP, 
albeit a few semi-hard particles, with T p E, still exist. For even larger times, these semi-hard 
particles will themselves disappear via democratic branchings, so there will be an increasing fraction of 
the total energy which is carried by the soft gluons with p < T. This trend is indeed visible in Figs. 9 
(f)-(h). However, one should not conclude that all this energy has already thermalized: the soft gluons 
which propagate together with their (semi-)hard sources along the light-cone z = t cannot be thermal. 
This is already illustrated by the last plot in Fig. 10: the momentum distribution corresponding to 
z = t = 2 trei is peaked at small p ^ T, yet it strongly deviates from a thermal distribution. 

To better distinguish between thermal and non-thermal (soft) gluons at late times, we have ex¬ 
hibited in Fig. 11 the z-distribution of the energy and number densities, defined as 

£it,z) = j dp\p\f{t,z,p), n{t,z) = jdpf{t,z,p). (5.5) 

As visible in these plots, even for times as large ast = 1.5tbr(F'), where we know (say, from Figs. 9 (f) 
and (g)) that the energy is preponderantly carried by soft quanta with p ~ T, the energy distribution 
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E=90T 



Figure 11. The energy density and the gluon number density are shown as functions of z at different times. 
The grey vertical lines indicate the location of the light-cone, that is, z = t, for each value of t. 


t — 141|.0|, z - 5.0 t|.@| 



Figure 12. The gluon distribution f{t, z,p) for E = 25T is plotted as a function of p ioi t = 14.2 t^ei — 3thr{E) 
and z = thriE) = 5trei and compared to a thermal distribution with the same normalization at p = 0. 

is still strongly peaked at z = t, meaning that most of these soft gluons are not thermal: they have 
been emitted at late stages and did not have the time to thermalize. This finding is in agreement with 
the discussion following Eq. (4.29), where we noticed that the flux of soft gluons is largest towards 
the late stages of the cascade. The situation changes at the larger time t = 3thr{E), that we can here 
access only for the jet with E = 25 T (by looking at t = 14.2 t^ei — 3tbr(25)). In that case, we see 
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that both densities, £(t,z) and n{t,z), peak well behind the light-cone — in particnlar, the number 
distribution peaks around z ~ thri^) = Streb agreement with Eq. (4.33). This strongly indicates 
that the entire glnon distribution produced by the jet has thermalized hy t = the jet is fully 

quenched. This conclusion can be also checked by plotting the distribntion f{t,z,p) as a function of 
p for t = 3tbr(-E) and, say, z = thriE): this is shown in Fig. 12 which indeed features an almost 
perfect thermal distribution. The distribution of such a fully quenched jet in longitudinal phase-space 
is illustrated by Figs. 8 (h) and 9 (h). Clearly, this is very similar to the late-time distributions found 
in Sect. 4, cf. Figs. 5 and 3: a distribution symmetric in 2 which extends via diffusion. 

5.3.2 Energy loss towards the medium 

Given our numerical results, as presented in the previous subsection, it is furthermore interesting to use 
them to extract the energy lost by the jet towards the medium. A priori, this involves two components: 
the energy dissipated into the medinm via the drag force (physically, this is the energy transferred to 
the plasma constituents through elastic collisions) and the energy taken away by the gluons from the 
jet which have reached a thermal distribution in momentum (since such gluons cannot be distinguished 
from the medium constituents anymore). As noticed at the end of Sect. 4.2.2, the drag component 
can be minimized by choosing p* = T, so it shonld be enough to compute the energy carried by the 
thermalized part of the gluon distribution. Still, to avoid any uncertainty concerning the contribntion 
of the drag, it is preferable to compnte the energy loss as the difference between the original energy E 
of the LP and the energy carried by the jet constituents which have not thermalized. This definition 
too is a bit ambiguous though, because the distinction between thermal and non-thermal gluons is 
not really sharp, as already noticed. Yet, we have seen that the gluons in the tail of the distribntion 
^ ^ ~ ^rei are approximately thermal, whereas those which belong to the front (z > t — Aei) are 
still far away from thermal eqnilibrium — at least for not too late time, t < ffj.{E), when the front 
still exist (see e.g. Fig. 10). This observation motivates the following definition for the energy loss via 
thermalization: 

roo POO 

A£’ther(t) =E- dz dppf{t,z,p). (5.6) 

Jpt 

For sufficiently large times t ^ ffriE), all the gluons \\e z < t — trei (the front disappears) and 
— E. But the most interesting situation in view of the phenomenology at the LHC, is that 
where the medium size L is small relative to the branching time tbr(-®) for the LP, hence AFther is 
small compared to E. 

In Fig. 13 we present our numerical results for AFther(t) as a function of t (or, eqnivalently, the 
medium size) for the two energies of interest, F = 25T and F = 90 T. For comparison, we also 
show the corresponding prediction AFfiow(i) of Eq. (4.29); this would be the energy transferred to 
the medium in the ideal case where the plasma acts as a perfect absorber for the gluons with p 
(without affecting the branching dynamics at p > T). Not surprisingly, the energy loss AFther(i) for 
the ‘physical’ cascade remains significantly lower than the ‘ideal’ expectation AFflow(i) for all times 
t < that is, so long as the gluon cascade has not fully thermalized. (For large times t 3> tbr(-£'), 

these two quantities approach to each other, as they both converge towards the total energy F, as they 
should.) The main reason for this discrepancy is the fact that the medium is not a perfect sink: there 
is a delay in the thermalization of the soft gluons produced via branchings and, as a result, gluons with 
p ~ T can still propagate at the speed of light over a time interval At > after their production. 
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Figure 13. The energy loss towards the medium AEtherit) (in units of T), Eq. (5.6), plotted as a function t 
for E = 25 T and E = 90T. The result is compared to the flow energy (4.29) which applies to an ideal cascade. 




Figure 14. The gluon distribution in momentum at t = z for two energies, E = 25 T and E = 90T, and for 4 
values of time, which are now chosen to be the same in units of tbr(^) for both energies. The figures show a 
rather broad window of approximate scaling behavior, / oc 1/p^/^, at not too large times t < t\^^{E). 


Since the production rate increases with time, cf. Eq. (4.29), it is natural that the difference between 
A£'flow(^) and AEther(^) increases as well so long as the LP still exists, i.e. for t < tbr(-®)- This trend 
is indeed visible in Fig. 13. 

This being said, the energy loss AEther that we have numerically found is significantly large. By 
inspection of Fig. 13, we see that A£'ther(T) ~ 30T (= 15GeV) for a jet with F = 90T (= 45 GeV) 
and for a medium size L = SFei (= 5fm). Furthermore, as also visible in Fig. 13, the energy loss rises 
quite fast with time and hence with the medium size L : at small times t <C tbrC-F), one roughly has 
AFther(i) oc in agreement with Eq. (2.6), whereas for larger times t > tbr(-F), AFther(i) approaches 
the total energy E of the LP: the jet is ‘fully quenched’. 

Whereas the results in Fig. 13 point out towards a failure of the hypothesis of a ‘perfect sink’, 
this failure remains quite mild, especially at early times t <C tbr(F'). This is already suggested by the 
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qualitative similarity between the gluon distribution produced by a ‘real’ jet, as shown in Fig. 8, and 
that generated by an ideal gluon cascade, cf. Fig. 3 and Fig. 6. To have a more quantitative test 
in that sense, we have studied the gluon spectrum near the front of the jet, that is, the distribution 
f{t, z,p) produced by the kinetic equation at z = t. At small times t < 0.5thr{E), the numerical results 
in Fig. 14 exhibit a relatively wide window at T < p E where the front distribution function shows 
the same scaling behavior, f{z = t,p) oc as the ideal branching process (compare to the curved 

‘w/o cutoff’ in Fig. 7). As repeatedly stressed, this scaling law is a hallmark of wave turbulence. Fig. 
14 should be contrasted to the corresponding results for the spectrum f{t,p) = f dzf(t,z,p) (the 
curves denoted as ‘full’ in Fig. 7), which show no scaling window at all. So, in this respect at least, 
the gluon spectrum is potentially misleading, as anticipated in Sect. 5.2. 

More generally, the ensemble of studies that we have performed in this paper demonstrate that 
the detailed phase-space distribution is better suited than the gluon spectrum for understanding the 
in-medium evolution of the jets. 

6 Conclusions and perspectives 

In this paper, we have presented a first study of the thermalization of the soft components of the 
gluon cascades generated via multiple branchings by an energetic parton which propagates through 
a weakly-coupled quark-gluon plasma. Our overall picture is rather simple and physically motivated, 
and in our opinion it is also quite robust: indeed, this picture is almost an immediate consequence of 
the strong separation of scales between the characteristic time for the medium-induced branchings of 
hard gluons and, respectively, the relaxation time for the thermalization of soft gluons. In trying to 
establish this picture beyond parameter estimates, we met with several difficulties and subtle points, 
for which we proposed at least partial solutions. 

A major difficulty is the overall complexity of the problem, that we have tried to circumvent via 
suitable approximations, notably by carefully separating the gluons from the jet from those in the 
medium and by projecting the dynamics onto the one-dimensional, longitudinal, phase-space. These 
approximations are fully justihed for the sufficiently hard gluons in the cascades, with momenta p ^ T, 
which control the dynamics of multiple branchings. On the other hand, these approximations becomes 
less justified when moving to the softer gluons with p < T, where they are at most qualitatively 
right. (But we have tried to carefully argue that they correctly reproduce the relevant time scales 
to parametric accuracy.) A particularly subtle approximation refers to the branching dynamics near 
the lower end of the cascades, at p ~ T, where we expect the cascade to terminate, on physical 
grounds. In our calculations, we have simply cut off the branching process at a scale p* ~ T and 
found that, thanks to the smearing effect of the elastic collisions, the results are not very sensitive to 
the precise value of this cutoff. (We have indeed checked that our numerical results remain qualitatively 
and even semi-quantitatively unchanged when varying this cutoff by a factor of 2 around its central 
value.) But it would be of course important, both conceptually and phenomenologically, to have a 
dynamical implementation of this cutoff, which in turn requires a consistent treatment of the full gluon 
distribution at soft momenta, including the inherent non-linear effects. 

This discussion points towards the many ‘technical’ limitations in our approach, with potential 
physical consequences, which will be hopefully lifted by more detailed, future, analyses. In principle, 
the theoretical framework for such studies is well defined: this is set by the general kinetic equations 
alluded to in Sect. 3, that can be found in the literature [36, 43]. A main difficulty as compared to 
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previous numerical studies of such equations in the literature [39, 51, 52], is that fact that, for the jet 
problem at hand, one needs to explicitly deal with the strong spatial inhomogeneity introduced by the 
hard components of the jet. 

Another interesting direction of research refers to a better understanding of the implications of 
the present picture for the phenomenology of jets in heavy ion collisions at RHIC and the LHC. Our 
first estimates for the energy loss via thermalization, which are of course very raw and must be taken 
with a grain of salt, are quite encouraging in that sense. In our opinion, it makes sense to compare, 
at least qualitatively, the quantity ALlther introduced in Eq. (5.6) (the energy carried away by the 
thermalized gluons in the tail of the jet) with the energy imbalance at large angles, as measured in the 
context of the di-jet asymmetry. The detailed analyses of the corresponding data, notably by the CMS 
collaboration [20, 26], demonstrate that the energy imbalance is carried by an excess of soft hadrons 
{pt ^ 2 GeV) propagating at large angles. It looks natural to associate these soft hadrons with the 
‘thermalized gluons’ in our current set-up. If so, it is interesting to notice that our estimates for AE^her 
in Fig. 13 are in the ballpark of 10 to 20 GeV, a value which is not unreasonable for the phenomenology 
alluded to above. But of course further studies, to remove some of our theoretical uncertainties and to 
better defines experimental observables, are still needed before aiming at a detailed comparison with 
the phenomenology. 
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